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Determination of Niobium and Tantalum in Titanium-Base 


Alloys 


John L. Hague and Lawrence A. Machlan 


A procedure for the determination of niobium and tantalum in titanium-base alloys is 


described. 
oxidized with a minimum of nitrie acid. 
containing 


a strong quaternary amine anion-exchange resin (Dowex-—1). 


The sample is dissolved in a mixture of hydrochloric and hydrofluoric acid, and 
rhe resulting solution is transferred to a column 


Titanium and 


other alloying elements, except niobium and tantalum, are removed by elution with an 


ammonium 


chloride-hydrochlorie-hydrofluorie acid 
elution with an ammonium cbloride-hydrofluoric acid solution. 
separate eluate of slightly acid ammonium chloride-ammonium fluoride solution. 


Niobium is removed by 
Tantalum is collected in a 
Boric 


solution. 


acid is added to complex the fluoride, the earths acids precipitated with cupferron, ignited 


to the pentoxides, and weighed. 
1. Introduction 


One of the difficult problems of analytical chemis- 
try has been the separation of niobium and tantalum 
from each other and from titanium when large 
amounts of titanium are present. This combination 
has become of practical importance with the develop- 
ment of MST grade 821 titanium-base alloy contain- 
ing 2 percent of niobium, 1 percent of tantalum, and 
8 percent of aluminum. 

It has been shown [3] * that titanium, niobium, and 
tantalum can be easily separated by anion-exchange 
in hyvdrochloric-hydrofluoric acid combinations. 
Some data are given here for the elution behavior 
of tin, molybdenum, vanadium, and iron, used as 
alloving elements in the titanium metal industry, 
and for phosphorus and tungsten because of their 
presence in naturally occurring ores or in concen- 
trates. While a question has arisen [1] as to whether 
it is possible to separate this group of elements from 
niobium and tantalum, it should be pointed out that 
present knowledge of the subject indicates that a 
fairly high concentration of hydrofluoric acid is 
required to eliminate interfering elements and to 
retain the niobium and tantalum on the columns for 
subsequent separation. 

The titanium metal allovs are dissolved in a mix- 
ture of hydrochloric and hydrofluoric acid and oxi- 
dized with a minimum of nitric acid. The solution 
is transferred to an ion-exchange column containing 
the strong tetramethylamine anion-exchange resin, 
Dowex-1l. The titanium, and any other interfering 
elements, are removed by elution with an ammonium 
chloride-hydrochloric-hydrofluoric acid solution and 
discarded. The niobium is collected in an eluate 
containing ammonium chloride and hydrofluoric 
acid. The tantalum is collected in a separate portion 
by elution with a slightly acid ammonium chloride- 
ammonium fluoride mixture. Boric acid is added to 
these eluates to complex the hydrofluoric acid, and 
the earth acids are then precipitated with cupferron. 
The cupferron precipitates are ignited to the pent- 
oxides and weighed as Nb.O; and TasQs. 

The procedure requires 2 to 3 days to carry to 


*Figures in brackets indicate the literature references at the end of this 


paper 


completion, as about 12 hr are required for the ion- 
exchange separations, and some time is required for 
the solution of the sample and the cupferron precipi- 
tations. 


2. Apparatus and Resins 


Ton-E-xchange Columns. The columns are con- 
structed of polystyrene and are approximately 12 in. 
long and 1 in.i.d. A simple column can be prepared 
from a 12-in. length of mon heer tubing as follows: 
pe tee 
The bottom of the tube is closed by a waxed No. 5 
rubber stopper with a *.-in. hole. A 6-in. length of 
polystyrene tubing, *is-in. o. d. and Me-in. bore, is 
inserted into the hole flush with the upper surface 
of the stopper. Another 6-in. length of this tubing 
is attached to the smaller tube with a 2-in. length 
of Tygon R tubing, and the flow controlled by a 
hose-cock clamp on the Tygon tubing. 

If a number of analyses are to be made, it is 
convenient to arrange the columns so that a number 
can be operated with a minimum of attention [5]. 
Plastic columns suitable for such an assembly have 
been developed [6], and are equipped with ‘‘Dole- 
type” fittings of polystyrene. Inlet and outlet 
tubes are of polyethylene; flexible connections where 
necessary are made of Tygon R tubing. The flow 
of solutions is controlled by hose-cock clamps on 
these flexible connections. The connections should 
be made and carefully checked in order to avoid any 
possibilitv of leakage of the solutions containing 
hvdrofluoric acid. An assembly of this type is 
illustrated in figure 1. 

Resins. Dowex—1, 200- to 400-mesh, S- to 10- 
percent divinylbenzene cross linkage are used. 
Experience over a period of several years has shown 
that the mesh-size of these resins may vary consider- 
ably from lot to lot. To avoid difficulty, the resins 
as received are air-dried, and sieved through a 270- 
mesh sieve. Material retained on the 270-mesh 
sieve is used for other purposes. Most of the fines 
are removed from the fraction passing the 270-mesh 
sieve as follows: Prepare a suspension of the resin 
in diluted hydrochloric acid (1+19).’. The coarser 

! Diluted hydrochloric acid (1+19) denotes 1 volume of cencentrated hydro- 


chlorie acid, sp gr 1.18, diluted with 19 volumes of water. If no dilution is speci- 
fied, the concentrated analytical reagent is meant 
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fraction is allowed to settle 10 to 15 min, and the 
fines removed by decantation. Repeat the process 
several times until most of the very fine material is 
removed from the suspension. Cover the bottom 
of the ion-exchange column with a ‘g- LO "s-IN. layer 
of acid-resistant vinyl chloride plastic wool. Add 
portions of the resin suspension to obtain a settled 


column of the resin 6 to 7 in. high. The loaded 
column should be run through several cycles of 
elution with alternate additions of diluted hydro- 
chloric acid (83+1) and diluted hydrochloric acid 
(1+9) to remove the remainder of the fines. The 
column is then washed with diluted hydrochloric 
acid (1+3) and is ready for use. Resin columns 
prepared in this way have been used for several years; 
the only maintenance required is to empty and 
refill the column with the resin charge if the flow 
rate becomes excessively slow due to packing. The 
resin charge should not be allowed to become dry. 

Polyethylene Ware. 250-ml and 600-ml Griffin 
form beakers are required. Bottles are used for 
preparing and dispensing acid mixtures containing 
hydrofluoric acid. 


3. Reagents 


Hydrochli ric-Hydrofluorie Acid Solution. Add 250 
ml of hydrochloric acid to 300 ml of water, then add 
200 ml of hydrofluoric acid, dilute to 1 liter with 
water, and mix well. This solution is referred to 
as 25-20 mixture in the procedure. 

Ammonium Chloride Solution (240 g/liter). Dis- 
solve 240 ¢ of ammonium chloride in 800 ml of 
water by warming, and dilute to 1 liter with water. 
Filter to remove insoluble material. This solution 
is used as a stock solution in preparing the three 
solutions which follow 

Ammonium Chloride - Hydrochloric - Hydroflucrie 
Acid Solution. Transfer 300 ml of ammonium 
chloride solution, 200 ml of hydrofluoric acid, 
and 125 ml of hydrochloric acid to a polvethylene 
bottle. Dilute to 1 liter with water, and mix well. 
This solution is referred to as 7-12 4-20 mixture 
in the procedure. 

Ammonium Chloride-Hydroflucric Acid Solution. 
Transfer 600 ml of ammonium chloride solution 
and 40 ml of hydrofluoric acid to a polyethylene 
bottle, dilute to 1 liter with water, and mix well. 
This solution is referred to as 14—4 acid mixture in 
the procedure. 

Ammon hum ('hloi ude -~Ammon hum bluoride Solu- 
tion. Transfer 600 ml of ammonium chloride solu- 
tion and 40 ml of hydrofluoric acid to a polvethvlene 
beaker. Adjust the solution to a pH of 5 to 6 with 
ammonium hydroxide (approximately 80 to 85 ml 
will be required), and dilute to 1 liter with water. 
This solution must be prepared with reasonable 
care, as a solution which is too acid will not com- 
pletely elute the tantalum in the volume specified 
in the procedure. A solution which is too alkaline 
will precipitate tantalum in the column, spoiling 
the determination being run and the one which 
follows. This solution is referred to as 14—4 neutral 
mixture in the procedure. 

(upferron Solution (60 q liter). Dissolve 6 ¢ of 
cupferron in 90 ml of water, dilute to 100 ml, and 
filter through a dry filter. This solution should be 
cooled (5° C) and prepared as needed, as the solution 


is not stable. 
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Cupferron Wash Solution. Add 25 ml of eupferron 
solution to 975 ml of diluted hydrochloric acid 
(1+9). Prepare as needed. 


4. Procedure 


As the samples dissolve quite rapidly, the rate of 
solution can be conveniently controlled by adding 
the sample to the cooled acid mixture. Transfer 
100 ml of the 25-20 mixture to a covered 250-ml 
polyethylene beaker. Add 5.00 g of the sample,’ 
in small portions, to the acid until the sample is 
completely transferred and dissolved. Warm slightly 
(40° to 50° ©) on the steam bath, and oxidize the 
titanium by the dropwise addition of nitric acid. 
Approximately 1 to 2 ml will be required, and the 
amount should be kept to a minimum because of 
the strong replacing power of nitric acid for niobium 
on the exchange column. Digest on the steam 
bath for 10 to 15 min, to remove the oxides of nitro- 
gen. Rinse the pl: istic cover and walls of the beaker 
with the 25-20 mixture, and dilute the solution to a 
volume of 150 ml with the same mixture. 

Transfer 50 ml of the 25-20 mixture to the column 
in small increments (5 to 10 ml), and drain the acid 
to 1 em above the resin bed, catching the eluate in a 
600-ml polyethylene beaker. Transfer the solution 
of the sample, prepared according to the directions 
in the previous paragraph, in small increments (5 to 
10 ml) to the column. Add the sample solution as 
it moves down the column until all the solution hes 
been transferred. Wash the beaker 4 to 5 times with 
5-ml portions of the 25-20 mixture, transferring the 
washings to the column. Finally wash down the 
sides of the column several times with approximately 
50 ml of the same solution. 

Add a total of 300 ml of the 7-12-20 mixture at 
a flow rate of approximately 100 to 125 ml/hr. Allow 
the solution to drain to the top of the resin, and wash 
the sides of the column with 4 or 5 portions (a total 
of about 25 ml) of the 14-4 acid mixture, allowing 
the solution to drain to the top of the column each 
time. Replace the beaker with another 600-ml! 
beaker, and discard the first eluate containing the 
titanium, vanadium, iron, tin, ete. 

Add a total of 300 ml of the 14—4 acid mixture, 
maintaining the technique and flow rate described in 
the preceding paragraph. Wash the sides of the 
column with approximately 25 to 30 ml of the 14—4 
neutral mixture in portions of 5 ml. Remove the 
beaker cont: 2ining the second fraction, and reserve 
it for the determination of niobium. Replace the 
beaker with a 600-ml polyethylene beaker. 

Add a total of 350 ml of the 14—4 neutral mixture 
and reserve the solution for the determination of 
tantalum. The column is cleaned by the addition, 
in increments, of 50 ml of diluted hydrochloric acid 
(1+3), and is then ready for the next sample. 


? Ifa smaller weight sample is used, approximately 25 ml of acid mixture should 

' used for each gram of sample, and an extra 25-ml portion added to maintain 
the high fluoride concentration required for a satisfactory separation of niobium 

Phi s point in the procedure provides a convenient and satisfactory place to 


Stop overnight if the elution is not carried through as a continuous operation 
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To the second fraction which contains the niobium, 
add 15 g of boric acid, 70 ml of hydrochloric acid, 
and 80 ml of water. Warm on the steam bath (30° 
to 35° C) to dissolve the borie acid. Cool to 5° C, 
and add slowly with continuous stirring 30 to 50 ml 
of cupferron solution.t. Add paper pulp, stir well to 
distribute the pulp, and allow to stand 10 to 15 min. 
Filter through a double thickness of 9-cm_ close- 
textured filter paper fitted to a Biichner funnel and 
precoated with a little filter pulp. Transfer the 
prec oyy to the funnel, and wash well with 400 ml 
of cold (5° C) cupferron wash solution. Transfer the 
precipitate pis paper to a weighed 25-ml porcelain 
crucible, and ignite at a low temperature until carbon 
is gone. Finally ignite to constant weight at 1,000° 
C, and weigh as niobium pentoxide. 

To the third fraction, which contains the tantalum, 
add 9 ¢ of boric acid, 85 ml of hydrochloric acid, and 
SO ml of water. Warm to dissolve the boric acid, 
precipitate, filter, wash, and ignite the tantalum as 
described for niobium. Weigh the tantalum = as 
tantalum pentoxide. 


5. Elution of Elements with Hydrochloric- 
Hydrofluoric Acid Solutions 


Of the elements that precipitate with cupferron 
[10], molybdenum and tin have proved to be the 
most troublesome to remove from the column (with 
a moderate volume of eluting solution) ahead of the 
niobium and tantalum. Both have been used as 
alloving elements for titanium metal, so elution 
curves covering a range of hydrochloric and hydro- 
fluoric acid concentrations are of interest. The 
curves for titanium, niobium, and molybdenum have 
been presented previously [3], and the curves 
obtained for quadrivalent tin are shown in figure 2. 
Curves for phosphorus (as phosphate) are given 
because it is present in ores and concentrates, and 
for iron (ferric) because of interest in ferroallovs of 
niobium and tantalum. As has been shown before 
(7, S], the retention of niobium and tantalum from 
solutions containing relatively high concentrations 
of hydrofluoric acid is good, so satisfactory separa- 
tions can be obtained in solutions of this type. 

Experiments demonstrated that the volume co- 
ordinate of the curves could be compressed somewhat 
by substituting ammonium chloride for a part of 
the hydrochloric acid. A “clean-up” solution con- 
taining 7.2 percent by weight of ammonium chloride, 
12.5 percent by volume of hydrochloric acid, and 
20 percent by volume of hydrofluoric acid was 
selected. The behavior of a number of possible 
interfering elements in this medium is also shown in 
figure 2. The excellent summary by Kraus and 
Nelson [9] of the elution behavior of twenty elements 
in hvdrochloric-hydrofluoric acid mixtures contain- 
ing from 0.5 to LN concentrations (approximately 
2 to 4% by volume) of hydrofluoric acid would 
by extrapolation, antimony, bismuth, and 


suggest, 


Blank determinations usually amount to less than 0.5 mg when the reagents 
conform to ACS specifications, and will ordinarily counterbalance the few tenths 
of a milligram of earth acid lost in the precipitation 
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thallium as potential interferences. Titanium alloys 
do not normally contain these elements, and they 
were not detected by spectrographic examination of 
the niobium and tantalum oxides obtained in the 
analysis of the alloys listed in table 2. Some earth 
acid ores contain considerable percentages of anti- 
mony and bismuth, which are normally removed by 
sulfide precipitation before the earth acids are 
isolated for separation. 

The behavior of some of the elements during ah 
analysis by the recommended procedure is shown 
by the first line of elution curves in figure 2. Titani- 
um, molybdenum, and tin are removed well ahead 
of the niobium fraction. The “tailing” of the 
tantalum fraction has received notice [2], and pre- 
liminary experiments indicate the pH of the tantalum 
eluting solution to be an important factor 


TABLE 2. 


6. Results 


The results obtained by the recommended proce- 
dure on a series of synthetic mixtures approximating 
the composition of the alloys to be analyzed are 
given in table 1. These values require little comment, 
being within the usual limits for gravimetric determi- 
nations. The titanium in these mixtures, available as 
commercially pure titanium metal containing neg- 
ligible quantities of niobium and tantalum, was 
dissolved in a hydrochloric-hydrofluoric acid mixture 
exactly as were the samples. Suitable quantities of 
niobium and tantalum pentoxides were weighed in a 
platinum crucible and fused in sodium bisulfate. 
The cooled fusion was dissolved in the acid mixture 
containing the titanium, the crucible removed, and 
the solution oxidized with nitric acid. 

This procedure was followed in preparing the mix- 
tures because the high-purity metals were not readily 
available. Aluminum was not added because it 
forms a sparingly soluble compound with sodium in 
the presence of considerable amounts of fluoride. 
Photometric procedures have been described [4] for 
the determination of small quantities (20 mg or less) 
of niobium and tantalum, and are not included in this 
paper. 


TABLE | Results of determinations of niobium and tantalum 


hy the recommended procedure in various sunthetic mirtures 


Niobium Pantalum 

Experi rita 
ment nin 

vided Added Found Ditfer Added = Found Ditfer 

ene ere 

mg mag ma¢ mg mq 

l 24. 5 24.4 .{ 175.0 174.7 0.2 

2 23.0 23.7 | 173.7 173. 5 ly 

3 " 135. 1 135. 1 0 20. 6 0. 6 0 

‘ 5 135. 2 135. 1 i 20.3 a4 X 9 

WS. WS. 4 2 WI 10. & l 

ti ; “3. 4 43.4 +0) ix Ss 18.7 ] 


The results obtained on one commercial and two 
experimental alloys are given in table 2. The de- 
terminations for each sample were obtained in two or 
three sets, and show satisfactory precision for gravi- 
metric determinations of these mixtures. One set of 
niobium and tantalum pentoxide precipitates for 
each sample was examined spectrographically for 
impurities. Moderate amounts (estimated at 0.1 
percent or less) of boron, silicon, and iron were usu- 
ally observed. The boron and silicon can be volatil- 
ized by treatment with sulfuric and hydrofluoric 


Re sults of determinations of niobium and tantalum hy the recommended proced cre in titanium-base alloys 


Niobium Tantalum 
Sampkh Sample 
nun. ner weight Nominal composition 
A veruet Runge No. of deter Average Range No. of deter 
value mination ilu minations 
WA-74 , 2.01 2. Wh te 2.02 md 0.0 te 1.00 - AL 8-Nb 2-Ta l 
WA-7 2.72 2.71; to 2.72 } 0). 34 $8; to C. 30 . Al 4-Nb 3-Ta 0.4-Sn 1 
WA-76 5 0. 58 0. 58» to 0. 5S { 3. 5S 3. 5Se to 3. 55 { Al 7-T'a 3.5-Nb 0.5-Mo 2-V 1.3 
7 } 
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acids, but experience over a period of several vears | 


has shown that this correction usually amounts to 
less than 0.2 mg. Iron is present due to impurities in 
the reagents, and is corrected for in the blank de- 
termination, amounting usually to no more than 0.2 
to 0.3 mg. Aluminum, calcium, magnesium, tin, and 
titanium were occasionally observed in amounts 
estimated to approach 0.1 percent, but usually in 
amounts of less than 0.01 percent. Occasionally, 
traces (estimated at 0.01 percent or less) of silver, 
chromium, manganese, and lead were reported. 
These impurity surveys, along with the data in table 
1 and figure 2, confirm the fact that the ion-exchange 
procedure furnishes a very satisfactory analytic “al 
separation. 


The authors are indebted to J. E. Griffin, Mallory- 
Sharon Titanium Corporation, Niles, Ohio; to M. J 
Miles of Titanium Metals Corporation, Henderson, 
Nev.; and to Galen Porter of Metals Research Labs., 
Electro Metallurgical Corporation, Niagara Falls, 
N. Y., who furnished the ingots for preparation of 
the samples WA-74, WA-75, and WA-76.  Par- 
ticular thanks are due to Sam Vigo, Watertown Ar- 
senal Laboratories, Watertown, Mass., whose efforts 
led to the formation of a task force (with J. L. Hague 
as chairman) to study the problem as a part of the 
analytical activities of the Metallurgical Advisory 
Committee on Titanium. 
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Properties of Zinc Borosilicate Glasses’ 
Edgar H. Hamilton, Roy M. Waxler, and Joseph M. Nivert, Jr. 


The glass-forming region of the system ZnO-B.O;-SiO, was surveyed. 


Glasses were 


made with compositions within the following limits in mole percent: 50 to 62.5 ZnO, 20 to 


$3.5 BOs, and 0 to 20 SiO». 


They can be melted and poured below 1,400° C. 
have low coefficients of thermal expansion, high values of elastic moduli and 


The glasses 
Poisson’s 


ratio, and refractive indices from 1.6409 to 1.6798. 


1. Introduction 


In a search for heat-resistant glasses, for which a 
low coefficient of thermal expansion is desirable, it 
was found that Winklemann [1]? reported a low, 
expansion zine borate glass containing 41 percent of 
BO, and 59 percent of ZnO by weight. He reported 
a coefficient of cubic expansion of 110107 over the 
temperature range of 10.35 to 92.88° C, which its 
equivalent to a coefficient of linear expansion of 
3.7*10-°. The glass has been reported to have 
other extreme properties including great hardness, 
large values of elastic moduli, and Poisson’s ratio 
(1, p. 193]. A glass-forming system, with a glass that 
is so highly recommended, seemed worthy of further 
study to determine if there were other glasses with 
similar or more desirable properties. 

K. Ingerson, G. W. Morey, and O. F. Tuttle [2] 
surveved the ZnO-B.O, and ZnO-B,O;-SiO, svstems 
and identified two zine borate compounds, ZnO-B,O, 
and 5ZnO-2B.0,, and a very large region of two- 
liquid immiscibilitv. They found that ZnO-B,O, 
has an extremely flat primary field—so flat that the 
determined values of liquidus temperature within 
the field were all the same, 1,000° C, within the 
limits of error of their determination. The ex- 
tremely flat primary field indicates a high degree of 
dissociation of ZnO-B,O; at the liquidus. In a 
number of other multicomponent systems the regions 
where stable glasses are obtained are in and adjacent 
to flat primary fields [3]. The ZnO-SiO, system was 
studied by E. N. Bunting [4] who identified one 
compound, 2ZnO-SiO,- 


2. Experimental Procedure 


The glasses were made in 500-¢ melts from batch 
materials of sufficient purity to satisfy the require- 
ments for optical glass. The melts were made in 
platinum crucibles and. stirred with platinum—10 
percent-rhodium propeller-type stirrers. The details 
of the melting procedure are published elsewhere 
[5, 6]. The sag points (accurate to +7° C) [7] of the 
glasses were used to establish annealing temperatures 
from which the glasses were cooled at 2.5° C per 
hour to 350° C. Refractive indices were measured 
with the Grauer refractometer [8] on small cubes cut 
from the annealed glasses. With this apparatus the 
refractive indices can be determined to +2«10~° 
without additional processing of the glass specimen. 

i The work described in this report was sponsored by the Wright Air Develop- 


ment Center, United States Air Force. 
2 Figures in brackets indicate the literature references at the end of this paper, 
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The liquidus temperature of each glass was deter- 
mined by a temperature gradient method [9]. The 
probable error of a single determination is +4° C, 
but the absolute error in the present data may be 
greater. Comparison with the data of Ingersonet al., 
to be discussed later, indicates that the error is less 
than 30° C. The thermal expansions of selected 
glasses were determined to within +0.1>< 107° by ar 
interferometric method [10]. Chemical durability 
was determined by the interferometric method of 
Hubbard and Hamilton [11]. Elastic moduli were 
measured by the dynamic method described by 
Spinner [12]. The accuracy of a single determination 
is £0.4 percent for Young’s and Shear moduli, and 
+1 percent for Poisson’s ratio. 

The compositions listed in table 1 are calculated 
from the batch formulations. 

The ZnO and SiO, contents of four of the glasses 
were determined by chemical analysis. The cal- 
culated and analyzed compositions of these glasses 
are compared in table 2. The results indicate that 
the changes in composition during melting of the 
glasses usually were less than 0.5 mole percent, but 
for a few of the glasses they may have been as high 
as 1 mole percent. 


3. Results and Discussion 


3.1. Glass-Forming Region of the System ZnO-B,O,- 


SiO, 


The compositions of the melts are given in table | 
and are plotted in figure 1. Clear binary zine borate 
glasses were made containing from 50 to 60 mole 
percent of ZnO. This composition range coincides 
with that of the extremely flat primary field of 
ZnO-B,O; as determined by Ingerson et al. [2]. 
All of the glasses could be melted below 1,400° C, 
The melts were very fluid and easily fined. 


3.2. Liquidus Temperatures 


The liquidus temperatures of many of the glasses 
are given in table 1. There is a gradual change in 
the composition of the glasses at and above their 
liquidus temperatures. The length of time that a 
sample is held in this temperature range will have a 
bearing on the liquidus temperature obtained, Suc- 
cessively higher values were obtained on repeat 
determinations. The liquidus values given in the 
table are the highest temperatures at which crystals 
were observed after a 30-min holding period. The 
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values are about 30° C lower than those r 
Ingerson et al. [2 


In either of these inve 





‘ported by 


stigations., 


failure to obtain equilibrium or losses due to volati- 
lization could account for such differences. 

All the clear glasses, free from ery stal defects on 
the surface or in the body of the glass, had liquidus 


temperatures between 900° and 1,110° C. 


A melt of the composition ZnO-B,O,, No. 40, was 


opaque at pouring temperatures. After 
was discovered that the cast block 
thin opaque shell surrounding a mass 


transparent glass. The ZnO contents of 


cooling it 


consisted of a 


of clear, 
the trans- 


parent glass and the opaque shell were determined 


with the following results: transparent 


sé *) 
glass 52.0 


mole percent, and opaque shell 0.9 mole percent 
So it appears that the region of liquid immiscibilits 
is slightly larger than indicated by Ingerson et al. 
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@, One liquid, cooled to glass; X, one liquid, cooled to glass plus crystals; A, 


two immiscible liquids 


(2) and includes the compound composition ZnO-B.Os. 
The temperature at which the melt is fined and 
poured will determine whether a clear glass or two 
phases are obtained. For example, repeat melts of 
the No. 40 composition produced two immiscible 
liquids when the melts were fined and poured below 
1,200° C; and when melted and poured at 1,250° to 
1 300° C, clear glasses were obtained indicating that 
at these temperatures only one liquid was present. 
When the clear glasses were heat-treated at tempera- 
tures above 957° C tiny globules appeared in the 
body of the glass, probably because of the separa- 
tion of an immiscible liquid phase. 

Zine borosilicate glasses containing less than 55 
mole percent of ZnO were characterized by the 
presence of thin films on the top surfaces of the cast 
blocks. These thin films produced interference 
colors when viewed by reflected light. The inter- 
ference colors were visible immediately after the 
melt was poured. The thickness of the films varied 
across the surface of each glass as well as from glass 
to glass. Presumably the thin films result from the 
presence in the melt of two liquid phases which are 
immiscible at pouring temperatures. 


3.3. Refractive Indices 


The refractive indices of the annealed glasses are 
given in table 1 and are plotted against the SiO, con- 
tent of the glasses in figure 2. The effect upon the 
refractive indices of the introduction of SiO, into 
zine borate glasses varies with their ZnO contents. 
For glasses containing 57.5 mole percent or more of 
ZnO, the refractive indices of the glasses increase 
with increase in SiQ, content and the rate of in- 
crease is determined by the ZnO content of the 
glasses. For glasses containing 55 mole percent or 
less of ZnO, increasing the SiO, content decreases 
the refractive index. 
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Figure 2. Change in refractive index with change in composi- 
w tion along ZnO isopleths. 


lhe valne adjacent to a line or a point indicates the ZnO content in mole percent. 


3.4. Thermal Expansion 


Three of the zine borate glasses had coefficients of 
thermal expansion in the temperature range of 25° 
to 100° C which were essentially identical to that re- 
ported for Winklemann’s glass. For the glasses 
containing SiO, the coefficients were even lower and 
approached that of Corning’s Pyrex Brand glass 
No. 7740. The coefficients of the glasses increased 
with temperature, but at higher temperatures were 
still considerably lower than for most commercial 
glasses which range from 7.0 to 10.0 107°, 


3.5. Elastic Moduli 


Elastic moduli were determined on four of the 
glasses by a dynamic method [12]. The specimens 
were in the form of bars ground to approximately 
554 by 1% by ¥ in. in size. With bars of this size 
and shape the resonant frequencies for flexure and 
torsion were between 1,600 and 2,000 eps and 
3,400 and 4,000 eps, respectively. The glasses have 
Young’s moduli above 800 kilobars * and Poisson’s 
ratios of about 0.3, both of which are higher than 
the values reported by Spinner [12] for common 
optical glasses. Poisson’s ratio for most glasses is 
usually between 0.2 and 0.25. 


One kilobar is equivalent to 1X 10° dynes/em?, 











TABLE 3. 


Composition 


Glass 
No. 
SiO» B,O ZnO 2 
Vole © Vfole * Vlole * 
12 40) 60 Severe t 
th 5 5A do 
2s 10 35 55 do 


* One fringe is approximately equivalent to depth of attack of 0.29 u 
> Severe: Attack so great it was impossible to determine the surface 

¢ ND: Not detectable 

1 SC: In addition to the indicated surface 


3.6. Chemical Durability 


ilteration there was severe 


Using an interferometric method [11] the chemical 
durabilities of three of the glasses were determined 
in solutions buffered at pH values of 2, 4, 6, 8, 10, 
and 12. The results are given in table 3. The addi- 
tion of SiO, to the zine borate glasses did not improve 
the resistance to attack by acids but did increase the 
resistance to strong alkali. The results indicate that 
the glasses should be serviceable in the pH range 
of 6 to 10, inclusive. 

A polished sample of glass No. 28 was exposed to 
alternate cycles of wetting and drying. Kach cvele 
was of 30 min. duration. During the wetting cycle 
the sample was heated from 60° to 65° C while the 
humidity was maintained at the saturation point, 
then the temperature was lowered to 60° C and the 
relative humidity to 84 percent for the drying cycle. 
After 96 hr exposure there was no detectable change 
in the appearance, light transmission, or weight of 
the polished plate. 


3.7. Surface Devitrification 


During annealing, the two component ZnO-B,O, 
glasses always developed a dull appearance on their 
exposed surfaces. The amount of this surface 
devitrification varied with the time and temperature 
of heat treatment. The three component glasses, 
containing as little as 5 mole percent of SiO,, always 
maintained a brilliant transparent surface during 
annealing. The tarnishing of the ZnO—B,.O, glasses 
during heat treatment is probably one of the reasons 
why Winklemann’s glass has not been adopted for 
commercial use. This tarnishing may be caused by 
a loss of BO, from the surface of the glass at the 
annealing temperature. 


4. Summary 


The glass-forming area of the system ZnO-B,O,- 
S10, was surveyed. Clear glasses are obtained in or 
adjacent to the primary field of the compound 
ZnO-B.O;. This primary field is extremely flat, 
indicating a high degree of dissociation of the primary 
compound ai the liquidus. Such a condition appears 
to be favorable for the formation of stable multi- 
component glasses. 
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Chemical durability of selected glasses 


alteration in fringes, * at pH 
*xposure, 6 hr at 80° C) 


t 8 10 12 
0.1 0.2 e¢ND Severe 
! 0.1 Do 
2 2 eND 40.1 (SC 


iteration in terms of fringes. 


ittack at the air-liquid interface 


The glasses have moderately low thermal expan- 
sions and high values of elastic moduli and Poisson’s 


ratio. 
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Design and Performance of a Block-Type Osmometer 
Donald McIntyre, G. C. Doderer,! and James H. O’Mara 


A block-type osmometer has been designed and constructed utilizing several new 
features which enhance the precision of the osmotic-pressure measurement and allow simpler 
manipulative techniques. The details of the construction as well as a discussion of the 
design and performance of this osmometer are incluced in the paper. 


1. Introduction 


An improved block-tvype osmometer has been 
designed and constructed that is capable of making 
precise measurements of the osmotic pressures of 
polymer solutions. This osmometer was designed so 
as to make the measurements as free as possible of 
determinate errors without necessarily complicating 
the laboratory manipulative work. It incorporates 
many of the desirable features obtainable with 
other osmometer designs. In undertaking this 
work it was realized that the semipermeable mem- 
brane is the ultimate controlling factor in any 
osmotic-pressure measurement, but it was also 
known that a good osmometer is essential in order 
to obtain a reliable evaluation of membranes. This 
hew osmometer was also constructed so that com- 
parative data might be obtained to evaluate the 
performance of commercially available osmometers. 


2. Construction of Block Osmometer 


An osmometer must provide two compartments 
which are separated by a semipermeable membrane. 
When a solvent is placed in one compartment and a 
solution containing an impermeable solute is placed 
in the other, equilibrium between these two liquid 
phases can be established by applying a pressure 
to the solution which has a lower solvent chemical 
potential. This pressure is called the osmotic 
pressure, and can be created by the transfer of 
solvent through the membrane to cause a height 
difference of the solutions in the capillaries con- 
nected to the compartments, or by the application 
of external pressure through a manometer. The 
present osmometer can use either method of pressure 
measurement. In figure 1, a portion of the machin- 
ist’s drawing for the osmometer is reproduced, and 
in figure 2, the assembled and open views of the 
osmometer are shown. The solution and solvent 
are introduced into each half of the osmometer 
through the standpipe, located at C, which is 
subsequently closed by the needle valve, located at 
D. The osmotic pressure is then measured as the 
height difference in the two capillaries, connected to 
the blocks at E. The wafflelike blocks that clamp 
the membrane were made of 440C hardened stain- 
The adapter, nut, needle, and standpipe 
were made of 308 stainless steel. The sleeve that 
joins the uniform-bore capillary through a graded 
seal to the adapter was made of Kovar. The needle 


less steel. 
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FIGURE lL. 
A, Guide pins; B, bolt holes; C, standpipe connection; D, needle valve seat; 
E, capillary connection. 


Drawing of osmometer. 


valves were made of a soft stainless-steel needle 
bearing against the hard stainless-steel block. 

Stainless steel was used as the basic material in 
order to make the osmometer useful for aqueous 
solutions, and also to obtain the necessary hardness 
for polishing the faces of the osmometer. Initial 
plans were to use a more corrosion-resistant but 
hardenable stainless steel in the 400 series, but the 
desired stainless steel was not available at the 
time. As a consequence the corrosion resistance 
is not as great as had been hoped for, but proved 
to be adequate after the whole surface was well 
finished and free of machining marks. 

Two osmometers were constructed. The first 
osmometer was hardened in an “inert’’ atmosphere 
at the National Bureau of Standards. The second 
osmometer was hardened at the United States 
Naval Gun Factory in a salt bath with an addi- 
tional second tempering step to avoid possible stress 
from hardening. The hardness of both was about 
52 to 54 Rockwell C. 

After the hardening operations, the first osmometer 
was ground by using a magnetic chuck, and the 
second was ground with a regular chuck in order to 
avoid the induced magnetization which makes the 
Internal parts more difficult to clean. However, the 














and half of 


Fully assembled osmometer (left an 


right). 


Figure 2. 
osmomeler 


magnetic chucking enabled a better polished surface 
and seems preferable for the best final results. The 
surfaces were finally ground planar within approxi- 
mately one-half of a fringe of white light over the 
whole area. An unexpected problem arose during 
the grinding of the flat surfaces. A concave surface 
was produced in the center of the osmometer because 
the small pillars decreased the resistance to the 
grinder and also because the heat treating tended to 
transform the outer laver into a harder austenitic 
phase than the center. 


3. General Performance 


All performance tests were conducted with gel 
“never-dried” cellophane No. 600 membranes. 
These membranes were very impermeable to solvent 
and had “zero-point” corrections in this osmometer 
of always less than 0.05 mm, and usually 0.00+ 0.02 
mm. <At first the Kovar metal sleeve that was 
immersed in the water bath began to corrode badly 
by pitting. After a few unsuccessful remedies, the 
Kovar was nickel plated and then coated with resin 
at the glass seal, where the plating was least effective. 
This proved to be successful. Perhaps a better 
method of sealing the capillary to the osmometer 
would have been the one used by Krigbaum and 
Flory [{1],? who formed a foot on the capillary itself 
and eliminated a glass-to-metal seal. 

Initially, the first osmometer was tested without 
removing the scale, which was caused by oxidation 
during the hardening process, from the outside of 
the osmometer. However, bubbles immediately 
began to form within the osmometer after its immer- 


? Figures in brackets indicate the literature references at the end of this paper 
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sion in water. This condition was accompanied by 
corrosion which never did cease. Later the scale 
was removed and the osmometer reassembled, but it 
was two weeks before the gas from the corrosion 
finally diffused out of the blocks. The second 
osmometer, which had no oxidized scale, never 
showed these effects. 

The filling and emptying of the osmometer was 
easily done with a long hypodermic needle with a 
blunt end. This procedure avoided any bubbles. 
Each side of the osmometer required about 7 ml to 
fill the block, although more was required to fill the 
standpipe and capillary. About 13 ml was used for 
each side, and the excess liquid was removed from 
each standpipe as soon as the measurement was 
underway. The needle valves were tightened with 
a torque wrench to 18 to 20 Ib-in. The osmometer 
did not need to be removed from the bath during 
the filling and emptying operations, so that osmotic- 
pressure measurements at nonambient temperatures 
were easily made without any severe dilatometric 
problems. The temperature variation of the bath 
was kept to +0.002° C in order to eliminate the 
thermometric effect in the 0.5 mm- and 1.0 mm- 
diameter capillaries used with the osmometer. 

Measurements were based upon the difference in 
height between the vertical solution and solvent 
capillaries using cathetometers that gave an accuracy 
of +0.02 mm. During the measurements on solvent 
permeabilities it was found that the sum of the heights 
of the liquids in the two capillaries of the block 
osmometer immersed in the water bath slowly in- 
creased with time, at about the rate of 0.2 to 0.4 
mm/hr. This was so slow that the membrane dis- 
symmetry (or even osmotic pressures) did not change 
with time once equilibrium had been reached. It 
was found that the liquid levels decreased extremely 
fast (6.5 mm/hr) when the osmometer was imbedded 
in a thermostatted box filled with cotton. These 
observations pointed to the fact that there was 
permeation of the solvent through the edges of the 
membrane. When the tests in the cotton-filled box 
were modified so as to have cyclohexane vapor in the 
area surrounding the membrane, the leakage rate 
dropped to practically zero. Nevertheless, after 
several weeks of operation with the osmometer 
directly immersed in water without any additional 
sealing of the membrane edge, the osmometer was 
opened for an examination of the membrane; it 
showed no traces of water permeation. The mem- 
brane apparently is securely protected against bath 
fluid penetration by the tight seal and by the addi- 
tional head supplied by the solvent in the capillaries. 

In order to prevent permeation through the edges 
of the membrane, a special gasket containing a metal 
foil encircling the membrane around the outer edge 
was used. The results obtained were similar to 
those obtained when polyvinyl alcohol was coated 
on the outer edge of the optical flats after the osmom- 
eter was assembled. This was tested when the os- 
mometer was inserted into a xvlene bath as well as 
in a silicone oil bath, with identical results. All 
subsequent measurements have been conducted in 
a silicone oil thermostat. With these precautions 
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the decrease in the sum of the capillary heights is 
0.2 mm/hr when 0.5-mm capillaries are used with 


eyclohexane. Undoubtedly some of this decrease is 
due to evaporation of the cyclohexane. However, it 
does seem necessary to have a gasket or sealing of 
some sort on any block osmometer to prevent slow 
permeations through the edge of the membrane. 


4. Osmotic-Pressure Measurements 


Osmotic-pressure Measurements were made with 
two different types of osmometer. The special block 
osmometer described earlier was tested and then 
used to evaluate the operation of a commonly used 
osmometer of different construction [2]. The Stabin 
osmometers performed well when they were assem- 
bled with the same type of membrane (gel cellophane) 
as that used for the NBS reference block osmometer. 
The membrane dissymmetries after 2 or 3 days were 
the same as those of the NBS osmometer, that is, 
less than 0.05 mm of toluene or cyclohexane. The 
Stabin osmometers were easily assembled and worked 
well from the beginning. A mercury seal was used 
around the positioning rod. 

Figure 3 shows the results of osmotic-pressure 
measurements on a polystyrene fraction in toluene at 
30° C with two Stabin osmometers and the new 
block osmometer. In this figure, h is the difference 
in heights of toluene in the two capillaries in centi- 
meters while C is the concentration in grams per 
deciliter. Figure 4 shows the same fraction of 
polystyrene in cyclohexane determined at 40° © 
with the two different types of osmometers. Figure 
5 shows the results of the osmotic-pressure Measure- 
ments on a fraction of polystvrene in cyclohexane at 
different temperatures using different osmometers. 

Table 1 indicates the speed with which solvent 
permeates the membrane in the respective osmom- 
eters. The constant, /, is determined from the slope 
of the curve represented in eq (1), h being the 
difference in column levels in centimeters at time, /, 
in seconds. 


In } I+ (1) 


Both types of osmometers achieved equilibrium os- 


motic pressures in evclohexane within 5 to 7 days, 
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FIGure 3. 
fraction in toluene at 30° C 


Osmoltic-pressure measurements on a polystyrene 
using the different types of os- 


mome\lers. 


A. NBS block; @, Stabin; @, Stabin 
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FrGuRE 5. Osmotlic-pressure measurements on a_ polystyrene 

fraction in cyclohexane at different temperatures using differ- 

ent osmometers. 

35°, Stabin osmometer; 40°, Stabin osmometer; 45°, NBS osmometer, 

although the NBS block osmometer approached 95 
percent of the equilibrium value within 1 or 2 days 
and the Stabin osmometer in 2 to 3 days. The solu 
tion concentrations for each osmometer were de 
termined initially by weight and also by differential 
refractive index measurements in an interferometer. 
After a run, both the solvent and solution in the 
block osmometer were checked interterometrically 
for changes or errors. The refractive indices in the 
NBS block osmometer checked within 5.6 107° 
units, which for cyclohexane-polystyrene solutions 
would correspond to a change in concentration of 
0.0033 percent if permeation of the polystyrene mole- 
cules alone accounted for the refractive index change. 
However, the limiting factor in such analyses is not 
the interferometry but rather the sma!l changes that 
occur in solvents upon being handled or being in con 
tact with metal. These changes are due either to 
absorption of contaminants or actual chemical re- 
actions which alter the solvent. 


TABLE 1. Permeability constants for cyclohexane 
Osmometer k 
cm/sec 
Stabin 1 3.3 1075 
Stabin 2 3.6 * 1075 
NBS block 1 3.8 107-5 


5. Discussion 
5.1. Design 


Osmometer design can basically be resolved into 
(1) the manner of supporiing the semipermeabl: 
membrane, (2) the critical dimensions of the cell, 
which determine the membrane surface to cel! 
volume ratio. (3) the mechanism of filling and 
emptying the osmometer, and (4) the method of 
measuring the osmotic pressure. Before discussing 





the design features in the NBS block osmometer it 
would seem necessary to point out that the design 
considerations were based upon the empirical knowl- 
edge that the final limit of error resides in the ‘‘zero- 
point correction”, or “membrane dissymmetry”’ 
This factor, even with a well-conditioned membrane, 
may be unknown to as much as +0.05 mm height 
difference of an organic solvent like toluene or cyclo- 
hexane. In view of this uncertainty, the measure- 
ment of the osmotic pressure need only involve a 
resolution of +£0.02 mm of solvent. 

The membrane is obviously best supported when 
both sides are firmly held in place to avoid ‘‘balloon- 
ing’ or other relaxation effects. Meyer [3], Fuoss 
and Mead [4], Flory [5], Carter and Record [6], 
and Trementozzi [7] have successfully used milled 
metal blocks to support the membranes. Various 
types of milled designs have been used. In spite of 
a slight loss in membrane rigidity the present in- 
vestigators undertook the more costly oper: ation of 
milling small pillars rather than circular grooves in 
order to gain the definite advantage of easier bubble- 
free filling. Demonstrations with a completed 
block, i half an osmometer, sealed with a glass 
plate have visually shown the complete avoidance 
of bubbles in the filling. 

Since the membranes are supported on optical 
flats it seemed that this tight seal of the membrane 
would make additional gaskets on the outer pe- 
riphery unnecessary. Other workers |S! have stated 
the need for such a gasket, and the present investi- 
gation has also demonstrated the need. even though 
other osmometers have been reported to operate 


successfully without the use of an additional seal 
14,5,7]. 
This type of membrane support has proved 


effective in avoiding drastic relaxation effeets in 
the membrane. The low membrane dissvymmetries 
and the quick approach of the zero pot to a con- 
stant value have been encouraging proof of this 
point. In addition it appears that the membranes 
currently being used are highly relaxed before they 
are introduced into the osmometer the low 
zero-point correction is quickly established after 
2 or 3 days and several solvent The devel- 
opment of the more rigid inorganic or glass mem- 
branes would, of course, eliminate the need of such 
precision finishing of the blocks, but after many 
years these membranes are still not yet widely used 
nor have they been critically tested. 

The surface area of the membrane that ts directly 
exposed to solution is approximately 41 em? and the 
volume of solution is approximately 7 ml, so that 
the surface/volume ratio is about comparable to 
that of other block osmometers. The effect or 
role of adsorption has not been studied enough to 
allow any more scientific choice of surface/volume 
ratios other than not making them too large. 

The ease of filling and emptying, as well as the 
ability to make an accurate concentration anaylsis, 
has been a real justification for the attempt to im- 
prove the block osmometer. The NBS model has 
shown that high quality workmanship can produce 
an instrument that has a low leakage and volatili- 
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zation rate with organic liquids, even those having 
boiling points as low as that of cyclohexane. De- 
vices built into the capillaries to avoid solvent 
evaporation have not proved necessary, even in 
work with cyclohexane at 45° C. 

The block-type osmometer has the added advan- 
tage that the pressure necessary for equalizing the 
solvent and solution activities may be easily ap- 
plied externally as well as through the difference of 
the heights in the capillaries. It is planned later 
to extend the range of osmotic pressures by adding 
external pressure to the capillaries. 

The heavy metal blocks of the osmometer are 
very useful in keeping a constant-temperature 
profile on both sides of the membrane throughout 
the region where the osmotic pressure is being 
created. A small temperature differential can mask 
the difference in activity between solvent and solu- 
tion, which is precisely the quantity that one is 
trying to equalize by a pressure increase. In the 
nonblock osmometers the lagging due to the inter- 
posed solvent partially compensates for this lack of 
metal, but certainly the block osmometer is to be 
desired for critical work. 


5.2. Comparative Measurements 


As the data in figures 3, 4, and 5 indicate, osmotic- 
pressure measurements using polymer fractions with 
the two osmometers are indistinguishable within 
experimental error. This makes routine osmotic- 
pressure measurements valid with either osmometer. 
However, for unfractionated materials a block 
osmometer becomes extremely valuable because it 
facilitates the detection of solute permeation. 

As an instrument for testing membranes and their 
permeability to solvent and solute species, the block 
osmometer can be extremely valuable. In addition, 
this osmometer can function easily as a null osmom- 
eter for dynamic measurements at any temperature 
compatible with the membrane. 


The authors express their thanks to W. Gnueg of 
the NBS machine shop and E. Webb and 
workers in the NBS optical shop for their cooperation 
and expert construction work. The authors also 
express their appreciation to L. H. Fawcett and his 
group at the United States Naval Gun Factory for 


his co- 


their cooperation in hardening the steel for the 
osmometers. 
6. References 
[1] W. R. Krigbaum and J. Flory, J. Am. Chem. Soe. 75, 
1775 (1953). 
(2) J. V. Stabin and E. H. Immergut, J. Polymer Sci. 14, 209 
1954) 
[3] K. H. Meyer, E. Wolff, C. G. Boissonnas, Helv. Chim 
Acta. 23, 430 +% 
[4] R. M. Fuoss and D. J. Mead, J. Phys. Chem. 47, 59 (1943). 
[5] P. J. Flory, J. Am. C hem. Soc. 65. 372 (1943). 
[6] S. R. Carter and B. R. Record, J. Chem. Soc. 1939, 660. 
[7] Q. A. Trementozzi, Monsanto Chemical Co., Springfield, 
Mass. (private communication). 
[8] A. F. Sirianni, L. M. Wise, and R. L. MeIntosh, Can. 
J. Research |B] 25, 301 (1947). 
WaAsHINGTON, October 28, 1958. 








ing 
Je- 
Pit 

in 


An- 
the 
Lp- 

of 
ter 


ing 


are 
ure 
yut 
ing 
isk 


lu- 


the 
er- 
of 
be 


ic- 
ith 
hin 
He 
er. 
ek 


it 


e1r 
ck 
on, 
n- 
ure 


of 
cO- 
ion 
lso 
his 
for 
the 








Journal of Research of the National Bureau of Standards 


Vol. 62, No. 2, February 1959 


Research Paper 2932 


Heat of Formation of Titanium Tetrabromide 
Raymond A. Nelson,' Walter H. Johnson, and Edward J. Prosen 


The heat of formation of solid titanium tetrabromide has been measured by the reaction 


of titanium metal with gaseous bromine 
the heat of reaction was obtained: 


Ti(c) + 2Br.(g) 


AHf? (25° C) 678.16 


4.60 ki/mole | 


a calorimeter. The following value for 


TiBry(c) 


162.08 + 1.10 keal/mole). 


Combination of this value with the heat of vaporization of bromine yields the standard 


heat of formation of TiBrs(c) from Ti(e) and Br,(liq) as AHf® (25° C) 


mole (— 147.40+ 1.10 keal/mole). 


1. Introduction 


Prior to this investigation the existing value for 
the heat of formation of titanium tetrabromide was 
that estimated by Brewer, Bromley, Gilles, and 
Lofgren {1]?, who gave —170 +10 keal/mole for the 
reaction of titanium with gaseous bromine. Their 
calculations were based on an old value for the heat 
of formation of titanium tetrachloride [2], and upon 
disproportionation studies and estimated heats 
of solution. 

Rossini, Wagman, Evans, Levine, and Jaffe [3] 
give —155 keal/mole for the heat of formation of 
titanium tetrabromide; this value was obtained by 
combining the value estimated by Brewer et al., 
with the heat of vaporization of bromine. 

The present investigation was undertaken to 
provide direct experimental data for the heat of 
formation of titanium tetrabromide. 


2. Source and Purity of Materials 


A sample of sponge titanium, obtained from the 
U.S. Bureau of Mines, was used for some of the 
preliminary experiments in this investigation. No 
information was available regarding the purity of 
this material except that it had a Brinell hardness 
of 136. 

A sample of commercial titanium sheet was used 
for other preliminary experiments. This sample 
was obtained from the Titanium Metals Corporation 
of America, who supplied the following analysis: 
99.74 percent titanium, 0.019 percent carbon, 0.047 
percent nitrogen, 0.08 percent iron, and 0.17 percent 
oxygen. 

For the final measurements, the results of which 
are tabulated in this paper, we used a sample of 
titanium prepared at the Naval Research Laboratory 
by the thermal reduction of titanium tetraiodide. 
This sample was cut from the same material as that 
used for measurement of the heat of formation of the 
tetrachloride [4]; the purity of this sample was in 
excess of 99.99 percent. 


Present address: U.S. Treasury Department, Alcohol Tax Unit, Washington 


25. D. C. 
2 Figures in brackets indicate the literature references at the end of this paper. 
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-616.72 + 4.60kj 


The bromine was taken from reagent grade ma- 
terial. A portion of the bromine was analyzed and 
found to contain 0.22 percent of chlorine, but the 
presence of iodine or of organic bromides was not 
detected. The effect of the chlorine is discussed in 
another section of this paper. 

A mass spectrometric analysis of helium made by 
the Mass Spectrometry Section of the Atomic and 
Radiation Physics Division showed that the helium 
contained less than 0.01 percent of air. It was 
purified further by passing it over titanium metal at 
600° C, copper oxide at 600°C, and through absorbers 
containing Ascarite, anhydrous magnesium perchlo- 
rate, and phosphorus pentoxide. 


3. Apparatus 


It was intended that the same calorimetric reaction 
vessel would be used for the experiments on titanium 
tetrabromide that was used for the tetrachloride [4]. 
However, since an operating temperature of 500° C 
was required, this vessel, which was designed for 
operation only up to 350° C, was not suitable. It 
was decided, therefore, to use instead a vessel previ- 
ously used for measurement of the heats of decompo- 
sition of diborane and pentaborane [5]. The titanium 
sample was placed between glass-wool plugs near the 
lower end of the quartz tube and was heated by the 
lower segment of the heating coil. A mixture of 
bromine and helium vapors entered the vessel at the 
bottom, passed over the hot sample and out through 
the glass helix. 

The bromine was contained in an external bubbling 
vessel which was weighed before and after each 
experiment. Valves were provided in the gas train 
whereby purified helium could (a) be passed directly 
into the reaction vessel, (b) be saturated with bromine 
and passed into the reaction vessel, or (¢) bypass both 
the bromine bubbler and the reaction vessel. 

The exit gases were cooled to the temperature of 
the calorimeter in the glass helix; upon leaving the 
reaction vessel they passed through a trap cooled 
in liquid nitrogen, where the excess bromine and the 
traces of TiBry vapor were collected. 

The calorimeter, thermometric system, and the 
apparatus for measurement of electrical energy have 
been described previously [6). 








4. Procedure 


In order to obtain a sufficient rate of reaction it 
was necessary to increase the surface area of the 
titanium by cutting the sample (0.2 to 0.3 g) into 
small pieces; a carefully cleaned bolt-cutter was 
used for this purpose. The sample was placed in 
the center of the heated portion of the vessel between 
plugs of glass wool. The vessel was then assembled 
and placed in the calorimeter together with 4,618 g 
of water. Moisture and oxvgen were removed from 
the vessel by evacuation followed by flushing with 
purified helium. 

The calorimeter temperature observations were 
divided into three parts: A 50-min reaction period, 
and two 20-min rating periods immediately preceding 
and following the reaction period [7]. Calorimeter 


temperatures were observed at 2-min_ intervals 
during the rating periods and at 1-min intervals 
during the reaction period. 

A mixture of bromine vapor and helium was 


introduced into the vessel at the start of the reaction 
period by bubbling helium through the bromine 
container. Two minutes later the vessel was heated 
electrically to 500° C and held at this temperature 
until the desired temperature rise of the calorimeter 
(1.5° to 2.0° ©) was obtained. The electric current 
was then interrupted, the bromine container by- 
passed, and the flow of helium through the vessel 
continued for about 20 min. The helium flow was 
then stopped and an additional 10 min was allowed 
for thermal equilibrium to be established. 

The rate of flow of helium was measured by means 
of a calibrated capillary flowmeter; the total quantity 
of helium was calculated from observations of the 
rate and time. The total quantity of bromine used 
was determined by weighing the bromine containet 
before and after each experiment ; the excess bromine 
was taken the difference between the total 
quantity and the stoichiometric quantity required for 
bromination of the sample. The temperatures of 
the helium and bromine entering the calorimeter 
were assumed to be the same as that of the calo- 
rimeter jacket through which they passed. These 
data were required only for the small heat 
capacity correction, where an error of 10 percent 
would affect the results by only 0.01 percent. 

The quantity of reaction was determined from the 
mass of sample, since the reaction was complete in all 
cases. The contents of the trap were analyzed for 
titanium in order to correct the results for the heat 
of sublimation of a portion of the titanium tetra- 
bromide. The solid titanium tetrabromide was dis- 
solved in a dilute sulfuric acid solution and the 
amount of titanium determined gravimetrically as 
the dioxide [8]. This analysis (accurate to about 
0.2) provided a check on the determination of the 
quantity of reaction. 

The calorimetric system was calibrated by per- 
forming separate experiments in the same manner 
except that the titanium and bromine were omitted. 
In one series of experiments the quantity of electrical 
energy required to produce a certain temperature 


as 
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rise is determined; in a second series this temperature 
rise is duplicated but a substantial amount of the 
required energy is derived from the chemical reaction. 
In this way the energy associated with a given quan- 
tity of chemical reaction is compared directly with 
an equal quantity of electrical energy as measured in 
terms of calibrated standards. 


5. Results and Calculations 


The results of the electrical energy calibration ex- 
periments are given in table 1, where £ is the elec- 
trical energy added to the system, g, is a correction 
for the heat carried into the system by the helium, 
ARe is the corrected temperature rise [7], Ae, is a 
correction for the heat capacity of the glass wool, and 
FE, is the electrical energy equivalent of the standard 
calorimetric system. The system was the same in 
each case except for differences in the quantity of 
glass wool. 


Results of the 


TABLE 1. electrical calibration ex pe iments 


Experiment I ] AR Ae r 
number 
Ohm ohm ram 
l $2278 6 $1 O. 161117 3.¢ AWSSS. 2 
2 39253 5 3. 8 1¥45939 3.f 2WSB51. 1 
3 31246 0 $4 Lala s 1 2415.8 
t 20594 1 22 147643 S 1 200450. 4 
5 31288 2 2.45 156234 0 20281. 1 
SH222 3 3. 4 LSOS5T 0 2003060 4 
ri 32135. 9 26 LHOBS5 0 2UUBSS. 4 
x 35748. 8 oy 178405 0 NO806 3 
u 7322.8 2 4 136364 28TH Y 
10 27373. 4 =e 130031 0 QOWOSH0. 8 
\Nlean YOOSES. | 
Standard deviation of mean 15. 6 


The results of the calorimetric experiments are 
given in table 2, where Ae is the deviation of the 
actual svstem from the standard calibrated system, ¢ 
is the total energy absorbed by the calorimeter, g 
is a correction to 25° C for the reactants, products, 
and the helium, and q, is a correction for the quantity 
of titanium tetrabromide sublimed during the experi- 
ment. The quantity g is obtained by the relation- 
ship: 

q—ARe(E, + Ae). 
The heat of reaction is obtained by the relationship: 
E—q+, 
moles ‘Li 


AH 4s, 


For these calculations the density of titanium was 
taken as 4.50 g/em® from the data of Fast [9]. The 
heat capacities of titanium, bromine, and helium were 
taken as 6.01, 8.60, and 4.97 cal/deg mole respec- 
tively, from the Selected Values of Chemical Thermo- 
dynamic Properties by Rossini et al. [3]. The heat 
capacity of titanium tetrabromide was taken as 28.93 
cal/deg mole from the data of Hall and Blocher [10]. 
The heat of sublimation of titanium tetrabromide was 
calculated as 14.9 keal/mole from the vapor pressure 








TasBie 2. Results of the calorimetric experiments on the reaction 
of bromine with titanium 
Ex- 
peri- - nog 
ment ARc Ag q I de q ri AIM 25°C) 
num- 
ber 
Ohm jiohm j j j j Moles kj/ mole 
1 ©0.131011 98 26251.7  23322.1 1.3 0.2 0.0043582 671. 95 
2 146616 5.6  20878.0 | 26009.9 2.1 2.1 . 0048326 678. 33 
3 132007 6.7 26450.9 23271.5 2.1 2.1 0046499 683. 76 
4 128886 625.0 | 25825.3 23462.0 2.4 2.2 0035004 675. 09 
) 149331 5.6 290918.0 26602.1 2.5 41.8 0049338 671. 04 
H 166368 3.2 33335.4 30075. 1 3.9 1.4 (048537 671. 30 
7 2505 9.4 $1210.6 37587.9 | 5.6 2.0 0052823 685. 14 
s 197521 10.4 39579.0 36071.4 4.8 L9 0050958 687. 76 
Mean...-.. 678. 16 
Standard deviation of the mean +2 34 


data of Seki [11]. All calculations were based upon 
the 1956 table of International Atomic Weights [12] 
For conversion to the conventional thermochemical 
calorie the following relationship was used: 


1 cal=4.1840 j. 


The results given in table 2 correspond to the fol- 
lowing process: 


Ti(c) +2Br,(g) 


TiBry(e) 


AHf°(25° ©) 678.164+4.60 kj mole 


162.08 + 1.10 keal/mole, 


If this value is combined with the heat of vaporiza- 
tion of bromine at 25° C, 7.34 keal/mole [3], the 
following standard heat of formation is obtained: 


Tile) + 2Br, (liq) = TiBry(e) 


AHf°(25° C) 616.72 + 4.60 kj mole, 


147.40+ 1.10 keal/mole. 


The corresponding heats of formation obtained 
from the preliminary experiments on sponge titanium 
and commercial titanium sheet are 146.9 and 

147.8 keal/mole, respectively. 

The stated uncertainty intervals are estimated 
over-all uncertainties, obtained by combining twice 
the standard deviation of the mean of the calibration 
and reaction experiments with reasonable estimates 
of all other known sources of error. 


6. Discussion 


A very slight amount of insoluble residue remained 
in the vessel after solution of the tetrabromide. This 
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residue could not be weighed accurately because it 
could not be separated from the glass wool. Approx- 
imate measurements indicated that the residue 
weighed less than 1 mg and it was assumed to be 
titanium dioxide. This quantity of material would 
account for approximately 0.1 percent of the original 
sample and would affect the heat of formation by 
only 0.04 percent. 

The presence of 0.22 percent of chlorine in the 
bromine would be expected to introduce another error 
of approximately the same magnitude. However, no 
evidence was found for the presence of the tetra- 
chloride among the reaction products. In any event, 
these errors are well within the uncertainty assigned 
to the heat of formation. 

Since the completion of this work there have been 
two independent determinations of the heat of forma- 
tion of TiBr,. Gross, Hayman, and Levi [13] meas- 
ured the heat of reaction of titanium with liquid bro- 
mine and obtained —148.1+0.25 keal/mole for the 
standard heat of formation. Schlifer and Schmidtke 
[14] made an independent determination by the same 
general method and obtained —147.87+0.04 keal, 
mole for the heat of formation. Both of these values 
agree with the results obtained in this investigation 
within the assigned uncertainties. 
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Adsorption of Polyesters on Glass, Silica, and Alumina’ 


Robert R. Stromberg, Alan R. Quasius, Samuel D. Toner, and Midgett S. 


Parker 


The adsorption of several saturated, linear polyesters has been studied under a variety 


of conditions. 
ester and poly(methyl methacrylate). 
and alumina. 
tained a large internal pore structure. 


The study also included di-n-butyl sebacate as a model low molecular weight 
The adsorbents were type E glass powder, silica, 
The glass powder and silica were nonporous materials; the alumina con- 
The rate of adsorption of the polyesters on glass 
was very rapid and the adsorpt on isotherm was irreversible. 
solvent was favored over adsorption from a good solvent. 


Adsorption from a poor 
A large number of sites are 


available on the adsorbent for bond formation and the adsorbed polymer most probably 


extends into the solution in a coiled configuration. 
appear to be involved in an electrostatic bond between the polymer and the glass. 


Metallic ions, in addition to silicon, 
The 


Langmuir equation is obeyed at high-equilibrium concentrations, although there is non- 


conformity at the low concentrations. 


1. Introduction 


Although the mechanism of adhesion has re- 
ceived attention for a number of years and con- 
siderable information has been obtained, no theory 
has been developed that is sufficiently general to 
apply to all systems and that will predict the 
behavior of new systems. 

Adhesion between two materials is restricted to an 
extremely narrow zone which involves a thin portion 
of each bulk material. Inasmuch as the locus of 
the interaction between a solid and an adhesive 
is at the interface, it is apparent that additional 
knowledge of this interface is required to determine 
the basic factors involved in adhesion. 

One approach to the determination of interfacial 
properties is a study of the adsorption of polymers 
onto solid surfaces and a study of the configuration 
of the polymer at the interface. The adsorption 
of polymers on solid surfaces is different in many 
respects from the adsorption of small molecules. 
The adsorption site and the configuration of an 
adsorbed small moelcule is usually fixed within 
known limits. The/inature and the adsorbability 
of the end groups ahd the “monomer” unit reoc- 
curring along the chain, together with inter- and 
intramolecular interactions, present a much more 
complex arrangement for polymer molecules. 

The work reported here is principally concerned 
with the adsorption of saturated linear polyesters 
onto glass, silica, and alumina. The effects of 
solvent, chemical structure of the polymers, and 
temperature were studied. 


2. Materials 


The adsorbents used for this study were glass, 
silica, and alumina. 

The glass was the kind used in the manufacture 
of glass fibers and is known as “type E”’. To ob- 
tain large surface areas and proper intermixing of 
"1 Presented before the Division of Polymer Chemistry of the American 
Chemical Society at the San Francisco meeting, April 1958. 


2 This work was done, in part, under a project sponsored by the Bureau of 
Aeronautics, Department of the Navy. 
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the adsorbent with the polymer solution, the glass 
was prepared in powder form. Type E glass mar- 
bles were broken in a “diamond” mortar, and ground 
in a ball mill containing type E glass marbles as 
balls. The powder was sieved through a U.S. stand- 
ard sieve number 325 with holes of 44 uw. A large 
quantity of glass powder, prepared in this manner, 
was thoroughly mixed to obtain homogeneity and 
stored under vacuum in sealed-glass tubes to mini- 
mize changes in the surface of the glass. 

An electron micrograph of a replica of the surface 
of this adsorbent is shown in figure 1. The lined 
ridges and grooves in the surface are the result of 
glass fractures that occurred in the mortar where a 
sharp blow brings about the size reduction. The 
rough, granular portion was caused by the grinding 
action in the ball mill. The roughness caused by 
the grinding contributes to the surface area by a 
factor of approximately two. The powder was non- 
porous and had a surface area, determined by nitro- 
gen adsorption and use of the Brunauer, Emmett, 
Teller (B. E. T.) equation, of 0.68 m?/g. The chem- 
ical analysis of this glass is given in table 1. The 
other adsorbents are also listed in this table. 


TABLE 1. Adsorbents 
Material Surface Particle size Chemical analysis 
area 
m2/g “ wt &% 
Silica *. 175 to 200 | 0.015 to 0.620 SiOzs, 99.0 to 99.7; Free 
1120, 0.2 to 2.0. 
Type E glass > (0, 68 <44 SiOs, 55; 
CaO, 16.0; 
AloOs, 14.5; 
B2O3,, 9.5; 
MgO, 5.0. 
Alumina 210 80 to 200 mesh 


® Information obtained from manufacturer. 
»b For chemical analysis of type E glass see reference 11. 


The silica was a type “prepared in a hot gaseous 
environment (1,100° C) by the vapor-phase hydrol- 
ysis of a silicon compound.” Some of the properties 
reported by the manufacturer are listed in table 1. 
It was noncrystalline and had a surface area deter- 














FIGURE I. 


rhe smooth, lined surface is produced by fractures resulting from the “diamond”? mortaring 


the ball mill 


mined by nitrogen adsorption and electron micros- 
copy of about 190 m*/g. The agreement between 
the values obtained by the two techniques indicates 
that the material was nonporous. The large surface 
area arises, naturally, from the small particle size. 
Inasmuch as large surface areas are experimentally 
easier to handle and adsorption onto nonporous 
structures is easier to interpret, the silica in this 
form provided a desirable adsorbent. 

The alumina was a type used for chromatographic 
analysis, and had a B. E. T. surface area of 210 m?/g. 
As this alumina had both a large surface area and 
relatively large particles, it appears that an internal- 
pore structure was associated with the material. 

The polymers studied are listed in table 2. The 
polyesters were furnished by James Farr of the 
Thiokol Chemical Corporation. They were pre- 
pared by refluxing the glycol and dibasic acid or 
anhydride for approximately 40 hr at temperatures 
up to 250° C with benzene or xylene added to remove 
the water of condensation. Following completion 
of the polymerization the systems were evacuated to 
remove the benzene or xylene and traces of water. 
The poly(neopentyl succinate) and poly(pentameth- 
vlene succinate) were crystalline materials, and the 
poly(propylene adipate) and poly(ethylene phtha- 
late) were amorphous. 

The poly(methyl methacrylate) was obtained from 
the Rohm and Haas Company. It had an intrinsic 
viscosity of 2.66 deciliters per gram in ethylene 
dichloride at 25° © and a calculated viscosity aver- 
age molecular weight of 1.2>< 10°. 


| approximately 5.8 g. 
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Electron micrograph of surface of type E glass particle. 


Che rough granular portion is caused by the grinding action in 


Solvents that would permit analysis of the polymer 
concentration by infrared absorption were required. 
Also, a comparison of adsorption from good and poor 
solvents was desired. Toluene was selected as a very 
poor solvent and chloroform as a good solvent. The 
solvents were dried by refluxing the toluene over 
sodium for several hours and by shaking the chloro- 
form over calcium chloride. Both were freshly dis- 
tilled prior to use. 


TABLE 4 Polume rs 
Mn» Min 
Poly (neopenty! Poly (trimethy 
succinate 1, 400 lene adipate 1, 000 
Poly (propy len Poly (pentameth 
idipate 5, 400 viene succinat 1 500 
Poly (ethylen Poly (ethylene 
phthalate 2, 000) idipate) 3, 800 
Di-n-butyl 
sebacate-. 314.4 
Poly (methyl — 
methacrylat« 1.210% M, 


* Nin, 


Number average molecular weight 


I,, Viscosity average molecular weight 


3. Procedure 


The amount of material adsorbed was determined 
by measuring the change in concentration of the 
polymer after equilibration of the solution. The 
concentration changes in the polyester solutions were 
measured by changes in the carbonyl absorption at 
A double-beam infrared spec- 
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A cell 
length of 0.4 mm was used and the polymer solutions 
were compensated with solvent in the reference 


trophotometer was used for these analyses. 


beam. A separate calibration curve, based on at 
least six polymer concentrations, was obtained for 
each experiment. The sensitivity was such that a 
change in absorbance of approximately 0.2 was ob- 
tained for each concentration change of 1 mg/ml. 
An experiment usually employed sufficient adsorp- 
tion tubes to complete one isotherm. 

The poly(methyl methacrylate) solutions were 
analyzed for polymer concentration by a nephelo- 
metric technique [1]... The polymer was precipitated 
from solution as a colloidal suspension by the drop- 
wise addition, with constant stirring, of 10 ml of a 
butanol-methanol mixture (2:1) to 2 ml of polymer 
solution. The amount of light scattered at 90° by 
this suspension was measured with a light-scattering 
microphotometer. A rigid time schedule was fol- 
lowed for the entire procedure. This made possible 
reproducible results. During the measurement of 
the turbidity the stirring was discontinued. The 
concentration of polymer was determined from 
calibration curves. 

The adsorption was carried out for most experi- 
ments with 20 ml of polymer solution and either 10 g 
of glass powder, 5 g of alumina, or 1 g of silica. 
After adjusting the adsorbent and polymer solution 
to the required temperature separately, for most 
cases the solution was mixed with the adsorbent, the 
Pyrex tube rapidly heated with a torch, sealed, and 
the adsorbent-polymer solution mixture rotated at 
approximately 9 rpm in a constant-temperature 
bath. In some cases the glass powder had been 
heated to 350° C in vacuum for 24 hr prior to addi- 
tion of the polymer solution. These outgassed and 
heat-treated samples were mixed with polymer 
solutions, both without any contact with the air and 
after exposure to air. The solution-adsorbent mix- 
tures were usually rotated for approximately 20 hr. 
Following the adsorption period the tubes were 
centrifuged and the solution removed for analysis. 

For desorption studies appropriate quantities of 
solution were removed after adsorption and pure 
solvent added to the tubes which were then resealed 
and rotated at constant temperature for 2 to 3 days. 
The supernatant solutions were then analvzed for 
polymer concentration. 


4. Results and Discussions 
4.1. General Shape of Isotherms 


The adsorption isotherm of poly(methyl methac- 
rvlate) from toluene solution onto type E_ glass 
powder is shown in figure 2. The shape of the curve 
is similar to that obtained by Jenckel and Rumbach 
[1] for this polymer on aluminum powder, quartz, and 
glass wool. There is no appreciable change in the 
quantity of polymer adsorbed after an equilibrium 
concentration of approximately 0.1 mg/ml has been 
reached. This rapid attainment of a plateau at 


§ Figures in brackets indicate the literature references at the end of this paper 
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Figure 2. Adsorption of poly(methyl methacrylate) from 


toluene on type E glass at 30° C. 


low-equilibrium concentrations has also been re- 
ported by Koral, Ullman, and Eirich [2] for the 
adsorption of poly(vinyl acetate) on iron and tin 
powder. The adsorption sites are readily occupied 
by the polymer at low-solution concentrations. In 
the competition between the adsorbent and solvent 
for polymer, the adsorbent has a strong affinity for 
the polymer and the poly(methyl methacrylate 
molecules are held to the adsorbent until active sites 
are no longer available. In the case of poly(methy1 
methacrylate) the nature of the adsorption isotherm 
appears to be independent of the adsorbent, at least 
for surfaces so far investigated. 

The polvesters studied yielded adsorption iso- 
therms quite different from those of poly(methy! 
methacrylate). Isotherms obtained at 30° C for 
the adsorption from chloroform and toluene solu- 
tions of several polyesters onto type E glass powder 
are given in figure 3. The polyester isotherms do 
not level off into a plateau as was the case with 
poly(methyl methacrylate). The quantity of poly- 
mer adsorbed continues to increase at the higher- 
solution concentrations. This type of isotherm is 
in agreement with a statistical theory of polymer 
adsorption developed by Simha, Frisch, and Eirich 





28 
» 
20 ——— 
> 9 
o 
E 
Go os 
D ~ 
> 12 x =" ~ 
a 
4 2) S 
ad © ——r* —& 
e I ie —s 7 a — 
4 - — ——"o 
e geo 





EQUILIBRIUM CONCENTRATION, mg/m 


Figure 3. Adsorption of several polyesters on type E glass 
at 30” C. 
idsor ption from toluene: ©, Poly(trimethylene adipate D, poly (pentamethy! 


ene succinate); @, poly(propylene adipate). 
idsorption from chloroform: ®, Poly(pentamethylene succinate); @, poly 
ethylene phthalate); @, poly(ethylene adipate). 
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tion of polymer at 30° C, 


, Adsorption of poly(trimethylene adipate) from chloroform; @, adsorption 
of poly(propylene adipate) from toluene; @, desorption of poly(propylene adi- 
pate) with toluene. 


(3 to 6], which predicts an initial steep rise and a very 
slow attainment of the asymptotic limit. The 
theory was primarily developed for low-surface 
coverage. As saturation has not been attained, 
even at concentrations as high as 10 mg/ml, as 
shown in figure 4, some of the polyester isotherms 
at the lower equilibrium solution concentrations may 
be in this low-surface coverage range. In the case 
of poly(methyl methacrylate), that position of the 
isotherm which corresponds to little surface coverage 
exists only at very dilute solution concentrations, 
which are experimentally difficult to measure. 


4.2. Dependence on Solvent and Water 


The quantity of material adsorbed was very 
dependent on the solvent and about 2 to 4 times as 
much polymer was adsorbed from the extremely 
poor solvent, toluene, as from the good solvent, 
chloroform. This solvent dependence is, in part, 
the result of competition between solvent and 
adsorbent for the polymer and the greater adsorp- 
tion from the poor solvent is in agreement with the 
theoretical treatment [3 to 6]. The curves shown in 
figure 3, obtained from the same solvent are all 
similar in shape, although the spread among the 
curves obtained from the toluene solutions is con- 
siderably greater than those obtained from chloro- 
form solutions. 

Measurements for varying periods of time demon- 
strated that the adsorption of the polvesters onto 
the type E glass powder was very rapid, reaching its 


equilibrium value within 15 see, which was the 
shortest time interval studied. Poly vinvl acetate) 
has also been reported 12} to adsorb in short time 


intervals (within an hour) on nonporous surfaces. 
Diffusion of the bulky polymer mass does not appear 
to be a major factor in adsorption onto such surfaces. 
Adsorption in porous networks presents additional 
problems relating to the ability of the polymer 
molecule to enter small pores and the availability of 
large portions of the surface area to such molecules. 
This will be discussed later, 
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trom toluene 


Figure 5. Adsorption of poly(etnylene aaipate) 


on type E glass at 30° C. 


, Glass heat cleaned at 350° C in vacuum prior to adsorption; @, glass not 


pretreated prior to adsorption. 

No significant difference in the quantity of polymer 
adsorbed as a result of a pretreatment of the glass by 
heating in vacuum at 350° C for 24 hr was measured. 
The results of one such experiment are shown in 
figure 5 for the adsorption of poly(ethylene adipate) 
from toluene. Traces of water were undoubtedly 
still present at the glass surface and it appears that, 
above these trace quantities, water does not play a 
significant role in the adsorption of polvesters. 


4.3. Reversibility of Adsorption and Dependence on 
Temperature 


Although relatively high equilibrium concentra- 
tions were required for maximum adsorption of the 
polyesters, the adsorbed materials were tightly 
bound to the surface. This is shown by the results 
of desorption experiments on poly(propylene adipate ) 
in figure 4 and poly(neopentyl succinate) in figure 6. 
The irreversibility of adsorption is shown at several 
concentrations on the isotherm in figure 4. In the 
case of poly(neopentyl succinate) only those samples 
containing the maximum quantity of polymer ad- 
sorbed in that experiment were used for desorption 
No significant quantity of polymer was desorbed 
from the glass for either polyester. 

After adsorption of polv(neopentyl succinate) on- 
to glass from toluene solution, dilution with chloro- 








form brought about the desorption of an appreciable 
6 
es > —f- ———v 
4 oe @ a 
— ities 
> at » — 
< - 
— -" 
a o ’ 
faa P 
E LZ ‘ 
ro) x ¥ 
w 08 A 
a) 
x 
ra 
4 
RIUM CONCENTRATION, mg/m 
FIGURE 6. Adso plion and desor plion of — poly(neopent yl 


E 


glass. 


succinale 


on ly pe 
, Adsorption from toluene at 30° C.; @, desorption at 30° C with toluene 


, adsorption from toluene at S86.8° C; @, adsorption from chloroform at 30° C; 
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quantity of polymer. This desorption with the 
better solvent indicates that the bonds between the 
polymer and glass are more physical than chemical. 

The effect of temperature on the adsorption of 
poly(neopentyl succinate) is also shown in figure 6. 
There was no significant difference in the adsorption 
from chloroform solution in the temperature range 
30° to 50° C and only one curve is drawn through 
the points obtained at both temperatures. A 
negative change in adsorption with temperature 
was observed when toluene was used as the solvent. 
Both positive-temperature coefficients [2] and posi- 
tive and negative coefficients [1] have been reported. 
It is possible that the decrease in the amount of 
poly(neopentyl succinate) adsorbed at the higher 
temperature may be a result of increased solubility 
of the polymer in the poor solvent. 


4.4. Adsorption on Silica and Alumina 


Considerably more polymer was adsorbed per 
gram of silica or alumina than per gram of 
powder. The adsorption onto silica from chloroform 
solution for polv(pentamethylene succinate) and 
poly (neopentyl succinate) is shown in figure 7. The 
shape of the isotherm on silica is similar to that 
obtained on type E glass for polyesters. For the 
poly(neopentyl succinate) a plateau region had not 
been reached at concentrations as high as 15 mg/ml. 


glass 
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FiGgtReE 7. Adsorption on silica and alumina from chloroform 
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. Poly(neopenty! succinate) on silica; @, poly(pentamethylene succinate) 
on silica; @, poly (neopenty! succinate) on alumina 
On this curve the points shown for the equilibrium concentration 4 mg/ml 
were obtained over a time period of 10 min to 20 hr. Virtually no change in ad- 


sorption occurred during this interval 





Although the silica and alumina both had approx- 
imately the same surface area, about three times 
the amount of polymer is adsorbed on the silica as is 
adsorbed on the alumina. The general shape of the 
curve on alumina, as seen in figure 7, is different from 
that obtained for polvesters on the type E glass and 
silica, and similar to that obtained for poly(methyl 
methacrylate). The group of points at 4 mg/ml in 
figure 7 was obtained from a set of measurements 
in which the adsorption time was varied from 10 
min to 20 hr. The porous alumina probably con- 
tained areas that are accessible only through pores 
that were too small to permit access by large ‘poly mer 
molecules. If appreciable portions of the internal 
surface area were available to the large polymer 
molecules, the rate of adsorption should have been 
considerably slower than that indicated in figure 7. 

slow rate of adsorption into the porous network 
of charcoal has been reported by Hobden and 
Jellinek [7] and Claesson and Claesson [8]. It 
appears, therefore, that the internal surface area 
of the alumina was unavailable to the polymer. 

The quantity of polymer adsorbed per unit of 
surface area is shown in figure 8 for poly(neopenty! 
succinate) for the adsorbents studied. A much 
greater surface coverage was achieved on the glass 
than on the silica or alumina. The small quantity 
adsorbed on porous alumina is similar to the rela- 
tively low quantities of polymer adsorbed on charcoal 





relative to nonporous adsorbents [1]. It is also 
possible that the polymer with sufficient time 
(months) might diffuse into the pores. 
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FIGURE 8. Adsorption of poly(neopentyl succinate) on several 


adsorbents from chloroform at 30° C. 


™, Adsorption on type EF glass; @, adsorption on silica; @, adsorption on 


alumina, 


In addition to silica, type E glass contains caleium, 
aluminum, boron, and magnesium oxides. Inasmuch 
as considerably larger quantities of polymer were 
adsorbed on the glass relative to the other adsorbents 
it appears that the positively charged metallic ions, 
in addition to silicon, present in the surface of the 
glass have some advantages for bond formation, 
probably with the oxygen of the ester group. This 
bond is probably electrostatic in nature. If the 
results reported by Jenckel and Rumbach [1] are 
calculated as a function of surface area, a considera- 
bly larger quantity of poly(methyl methacrylate) was 
adsorbed on aluminum powder, which they report 
was most likely alumina on the surface, than on glass 
or silica. The amount adsorbed on glass was 
slightly more than on silica. These results can also 
be interpreted as indicating the importance of sur- 
face atoms such as aluminum for the formation of 
bonds with the adsorbed polymer. 
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Fiacture 9. Adsorption of di-n-butyl sebacate on silica from 


loluene at 30° C. 


To determine how a very low molecular weight 
material with structure similar to the polvesters 
would behave, the adsorption of di-n-butyl sebacate 
onto type E glass and silica was attempted. No 
adsorption was observed on the glass, using the 
techniques described earlier. An isotherm was ob- 
tained from toluene solution onto silica, however, 
and is given in figure 9. The shape of the curve is 
typical for the adsorption of polyesters, although the 
weight adsorbed was considerably less than that ob- 
tained with the polymers. No di-n-butyl sebacate 
was adsorbed from chloroform solution. 

Ellerstein and Ullman [9] reported very small or 
ho adsorption of a linear dimer on glass and attrib- 
uted this to the small entropy change, AS, that would 
occur upon adsorption of the dimer compared with 
the much larger change accompanying adsorption of 
a polymer. The adsorption of the small molecule 
would only release one solvent molecule into the 
solution while the adsorption of a polymer molecule 
on several sites would release several. If there is no 
appreciable difference in the heat of adsorption, A/7, 
for the dimer molecule compared to the polymer 
molecule, it is possible that this difference in AS 
accompanying the adsorption process could cause a 
positive free energy change, Af, for the dimer and 
a negative AF for the polymer. If the configuration 
of the di-n- -butvl sebacate in toluene is such that 
adsorption at more than one site occurs, permitting 
a larger entropy change, a negative free energy of 
adsorption may result, again assuming no change in 
All. A more open structure in chloroform may 
permit adsorption at only one site with no advan- 
tageous change in the number of degrees of freedom 
no appreciable adsorption. 

The availability of sufficient surface area tor the 
di-n-butyl sebacate molecule to lie flat can be de- 
termined from the isotherm for this material. If it 
is assumed (from a Fisher-Hirschfelder model) that 
the area occupied per molecule is 620 A’, then at the 
equilibrium concentration of 1 mg/ml, approximately 
1.1 10% were adsorbed and would 


and, consequently, 


molecules em 
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occupy about 70 percent of the available area if they 
were flat on the surface. At saturation, which is 
assumed to be approximately 35> 107° moles/m?, 
2.110" molecules/em*? were adsorbed and would 
occupy 130 percent of the surface, indicating that a 
portion of the molecule extended into the solution, 
at least at this concentration. 


4.5. Configuration of Adsorbed Molecules 


The deposition of adsorbed polymer in sequences 
separated by bridges was suggested by Jenckel and 
Rumbach [1] and an isotherm based on this concept 
has been derived by Simha, Frisch, and Eirich 
[3 to6]. In the case of di-n-butyi sebacate, loops may 
not occur, although sections of the molecule un- 
doubtedly extend into the solution in a manner 
similar to the beginning of bridges in polymers. If 
the molecule is adsorbed at only one location, some 
indication of the number of available sites can be 
determined, in this case sites sufficient for at least 
2.110" molecules/em?. 

The solvent and temperature dependence of the 
adsorption of polyesters can be explained in part by 
differences in solubility and in part by differences 
in the configuration of the polymer molecule. The 
decrease in free energy that occurs with adsorption is 
dependent on the solubility of the adsorbate. An 
increase in solubility as a result of either better 
solvent or increase in temperature would bring about 
a decrease in adsorption. In addition, the polymer 
chain probably exists as a coil with a radius of 
gvration that increases, within limits, with increased 
solubility. The increase in size of the polymer coil 
would also result in a decrease in adsorption. 

If the adsorbed polymers were to lie flat on the 
surface, the number of lavers can be calculated. 
The area of a “segment” of polv(neopentyl suc- 
cinate) is approximately 70 A*®. On the glass ad- 
sorbent there would be about 5 layers after ad- 
sorption from toluene and approximately 2.2 lavers 
from chloroform solution. There would be about a 
monolayer of this polymer adsorbed on silica and 
only 0.2 of a laver on alumina after adsorption 
from chloroform solution. 

Attempts to fit the data to the simple adsorption 
equations of Freundlich and Langmuir were made. 
The Freundlich equation was obeved over the 
vreater portion of the concentration range studied, 


but deviated at both high- and low-solution con- 
centrations. The Langmuir equation was obeved 


for poly(methyl methacrylate), except for the very 


low concentration range. There was marked de- 
viation at concentrations lower than 1 mg/ml for 


the polvesters. The Langmuir equation is derived 
for rigid molecules where the rate of adsorption is 
in equilibrium with the rate of desorption, and would 
not necessarily be expected to apply to polymer 
adsorption, especially the irreversible adsorption 
of poyesters to glass. Therefore, the deviation for 
the poly(methyl methacrylate) and the polyesters 
is not surprising. 

The theory developed by Simha, Frisch, and 
Kirich [8 to 6] predicts that a high surface coverage 

















where there is interference among the polymer 
segments, the derived isotherm approaches the 
Langmuir isotherm |2|. Both the polyesters and 
polv(methyvl methacrylate) appear to obey the 
Langmuir equation at the higher surface coverage. 
The polvesters depart from linearity at higher 
solution concentrations than the poly(methyl meth- 
acrylate). This may be a result of the more gradual 
approach to total-surface coverage for the polyesters 
where concentrations as high as 1 to 2 mg/ml are 
required before large-surface coverage and polymer 
segment interference is obtained. Koral, Ullman, 
and Eirich [2| found that the adsorption of poly 
(vinyl acetate) also appears to obey the Langmuir 
equation in the high-concentration range. Hobden 
and Jellinek [7] and Jellinek and Northey [10] 
reported that the adsorption isotherms of poly- 
styrene on charcoal obeved the Langmuir equation 
very well over the entire range they studied. 





5. Summary 


To summarize the significance of the adsorption 
data, it has been shown that the adsorption of 
polyesters onto glass is very rapid and that strong 
bonds, probably electrostatic in character, are 
formed between polyester and glass. This is demon- 
strated by the irreversibility of the adsorption and 
the larger quantities of polymer adsorbed on the 
glass compared to the silica. The adsorption from 
toluene, a poor solvent, is favored over adsorption 
from chloroform, a good solvent. This, together 
with the large number of sites available for bond 
formation, indicates that the configuration and the 
radius of gvration of the polymer molecule in solu- 
tion is an important consideration in determining 
the quantity of polymer that will occupy a given 
surface area. The number of lavers of polymer that 
would be adsorbed on the surface of the substrate 
would be of the order of 2 to 5 in this study if the 
molecule were to lie flat. As such an arrangement 
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would assume very large alterations in molecular 
configuration upon adsorption, it is reasonable to 
assume that the molecule must exist as a coil which 
extends out into the solution, contacting the surface 
at several points. This coiled molecule, therefore, 
also provides opportunity for intermolecular inter- 
action and entanglement during the adsorption 
process. 
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Wharton who conducted the majority of the infrared 
analyses. They also express their appreciation to 
Frank Reinhart, Sanford Newman, Bernard Ach- 
hammer, and to Robert Ullman (Polytechnic Insti- 
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Theory of the Effect of Drag on the Orbital 


Inclination of an Earth Satellite 
John P. Vinti 


The rotation of the earth’s atmosphere with the earth gives rise to a change in the 
inclination of the orbit of an earth satellite. In this paper it is assumed that the drag is in 
the direction of the air velocity relative to the satellite and that its magnitude diminishes 
so rapidly with altitude that it is appreciable only near perigee in an elliptic orbit. With- 
out further assumptions about the properties of the drag function, results are then deduced 
for the secular changes in inclination up to, but not including, the final ballistic stage of 
4 motion. 








1. Introduction 


lili ae 


Although the drag on an earth satellite is not well known, it is possible to express the 
secular changes that it produces in the inclination of the orbit, in terms of easily measurable 
quantities, provided only that the drag is given by a certain very general law. The present 
paper, which was suggested by a note by Wildhack,' gives a more complete treatment of the 
problem, with more liberal use of mathematical methods. 


2. Forces and Frame of Reference 


If we adopt a Cartesian reference system with origin at the center of the earth and axis 
directions fixed relative to the fixed stars, the following forces act on a satellite: 


F,— inverse square gravitational force of the earth, 
“> F, correction to F,, arising from oblateness of the earth, 


4 


F,=gravitational force of the sun, 
F ,, gravitational force of the moon; 
F,,= inertial force arising from acceleration of the earth, produced by the sun, 
F,,,, inertial force arising from acceleration of the earth, produced by the moon, 
F,=drag force. 

There are other forces, all of which may be disregarded as extremely small. Moreover, 
F.+F,, and F,,+F,,, are small compared to F,. We are thus left with the forces F,, F,, and 
F ,. 

The oblateness force F, has a secular effect on the orbit, inasmuch as it produces a secular 
precession of the line of nodes and of the line of apsides. In considering changes in the in- 
clination of the orbit, however, we are concerned with nutation and the effect of oblateness is 
only to produce a small oscillatory nutation. We shall therefore be concerned only with the 
forces F, and F,, up to that point where we find results dependent on the (oblateness-induced) 
motion of the perigee. At that point it will then be sufficient to quote results by Garfinkel? 
to find how oblateness influences the secular effect of drag on the inclination. 


3. Fundamental Theory 


With the above choice of reference system the center of the earth plays the role of a fixed 
point and the angular momentum p of the satellite is given by 


p mr» r. (1) 
Here m is the mass of the satellite, r is the radius vector to it from the center of the earth, 
and Fr is the velocity of the satellite. 


W. A. Wildhack, Science, 128, 309 (Aug. 8, 1958 
B. Garfinkel, Astron. J. 63, 88 (March 195s). 
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As an orthonormal vector basis we choose k pointing northward along the axis of rotation 
of the earth, i pointing in such a fixed direction that at time t=0 its extension pierces the 
Greenwich meridian, and j so that j=kX 1. 

The orbit is a space curve, but it is natural to define the instantaneous plane of the orbit 
as the plane of two successive radius vectors, i. e., as the plane of r and r or the plane perpen- 
dicular to p. It is then convenient to define the directed unit normal to the orbit as a unit 
vector 1, along p. The sense of the motion of the satellite relative to 1, is then given by the 
right-hand serew rule. We then define the inclination @ of the orbit as the angle between k 
and 1,. (In terms of observable quantities @ is thus the greatest latitude reached by the sat- 
ellite in a revolution.) Then 

p= pl,, (2) 
where p= p and 
1,=i sin 6 cos @+j sin 6 sin ¢+K cos @ (3) 


the angle @ being the angle between i and the projection of p on the equatorial plane i, j 
Then p, 6, and @ are the polar coordinates of p. 


The vector 


1y—k 1, ese 6 (4) 


is a unit vector lving both in the equatorial plane and in the plane of the orbit, so that it points 
along their intersection, the line of nodes. If 1, is a unit vector along r, it follows from (1) 
that 1, is perpendicular to 1,. Moreover ly is perpendicular to 1,, by (4), and so is 1, Ly. 
Thus the unit veetors 1,, ly, and 1, 1y are coplanar, since they are all perpendicular to 1,. 
If now y¥ is the angle between the radius vector and the line of nodes, so that 


y (ly, 1,), (5) 
we have 


1 ly cos y+1,&1y sin y, (6) 


where the positive sense for increase of y is from ly towards 1, Ly. 
Krom (3) and (4 
1, isin d+j cos @ (7) 


kx<1, i cos ¢—j sin ¢. (8) 
(‘omparison of (8) with (3) shows that 
1, k cos 6—k >» in sin @ (9) 


By (4) it follows that ly is perpendicular to k, so that k, iy, and kx 1, form another right- 
handed orthonormal basis, which turns out to be useful. 


4. Changes of Angular Momentum 
If F is the total force acting on the satellite, then 
p—r&F. (10) 
The secular changes with which we are here concerned are produced only by F,, so that for 
our present purposes 
p—rX<F,. (11) 


To find the rates of change of p, @, and ¢, it is now convenient to express each side of (11) in 
terms of the basis k, ly, and kx Ly. 
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The left side may be so expressed by means of (2), (3), (7), (8), and (9), with the result 
p=k(p cos @—p sin 0)6- lyp sin 66+kX1y(—p sin 6—p cos 60). (12) 
5. Drag Force and Drag Torque 
We assume the atmosphere of the earth to be rotating with it like a rigid body. If &, is the 
angular speed of rotation of the earth, the velocity w of the atmosphere (w for “wind’’), relative 
to the inertial system i, j, kK, is 


» 


w w,k xP. (133) 
The velocity of the satellite relative to the inertial system isr. Its velocity v, relative to 
the atmosphere is then 

vV,—r—w. (14) 


We assume the drag force to be in the direction of —v,, with magnitude dependent on r and 


the relative speed r—w. Thus 
F,=f(r,t—w )(w—r). (15) 


Equation (15) is the “very general law’? mentioned in the Introduction. If the drag is quad- 
ratic in the air speed r—w, then the function f is proportional to the air speed. No such 
special assumption is needed, however, for the purposes of this paper. 

The drag torque is now given by 


r<F,=frx<w—frxr. 16) 
Now by I 
r<r—p/m, (17) 
by (133 
r< w=w,r"*|k—I,(1,-k)], IS) 
and by 6) and (9) 
l.=KkK sin @ sin ¥+ly, cos ¥+kK& 1, cos @ sin y. (19) 


From (16) through (19) it then follows that 
rx<F, fp/m + feo.*{|k(1 —sin’6 sin? y)—1,y sin 6 sin ¥ cos Y—k™& Ly sin 6 cos 6 sin? YJ. (20) 
This is the expression for the drag torque. 

6. Basic Equations 


On inserting (12) and (20) into (10) and equating corresponding coefficients of the basis 
vectors kK, ly, and kx Ly, we find 


p cos 6—p sin 60 fm~'p cos 6+-fw.r?(1—sin? 6 sin? p) (21.1) 
p sin 6+ p cos 66 fm~'p sin 6+-fwe,-* sin @ cos @ sin*y (21.2) 
p sin 0 for? sin @ sin py cos y. (21.3) 
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Multiplication of (21.2) by <=. —1 and addition of the resulting equation to (21.1) gives, after 


division by expe: 





p ipo fm ‘p+fw.r*|cos* y exp( 10) +-cos 6 sin? y}, (22 
from which there follow 
p fm ‘Pp +-fw,r- cos 6 (23.1) 
po for cos” y sin @. (23.2) 
Equations (23), together with (21.3), are the basic equations for the changes in p, 6, and @ 
that arise from the drag. 
7. Stages of the Motion } 


We assume that the motion of a satellite takes place in three stages: 

Stage (1): The motion is that of a precessing and shrinking ellipse, with diminishing 
eccentricity, In which the elliptic elements are slowly varving functions of time. For our 
present purposes the relevant elliptic elements /, e, and WV ps along with the related quantities 
Pp, ’p, and r,, may be chosen as follows: 

p=magnitude of the orbital angular momentum 
/—semi-latus rectum 
é eccentricity 
r, radius at perigee 
} radius at apogee 
y,— Value of y at perigee. 
The elliptic orbit is then characterized by the equation 


l1+-e cos(¥—y,)]|-'. (24 + 


~ 


where 
/ (GAMm?) yp (25 
and 


p=mi “yp. (26 


M being the mass of the earth. From (24) it follows that 





/ (1 é) (2 
} /(] é) IN 
from which 
(=2r y.(f-+?o)' (29.1 
( / r,dQ / (29.2 


Stage (2): A quasi-steady spiral stage, in which there are no apsides, so that 7 is small 
and hegative throughout the stage, and in which the speed ris given approximately D\ the 
formula for a circular orbit without drag. With neglect of oblateness this relation is 


r GA\I rr. (30 
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Stage (3): A final “‘ballistic’”’ stage, which consists of a rapidly diminishing spiral, and 
in which v is not given approximately by the formula (30) for a circular orbit. Quantitative 
treatment of such a final stage would require an accurate knowledge of the drag function, so 
that we omit consideration of it in this paper, restricting our attention to the changes in 
inclination that may occur in the first two stages. 

We let 7,9 and 7,9 be the initial perigee and apogee radii in stage (1) and’, their final common 
value at the end of stage (1). The quantity 7, is then the initial radius in stage (2), which we 
take to have a final radius 7... In any actual case the quantities 7,5 and 7,9 would be given as 
initial conditions, while the radii 7; and r, would have to be determined by observation—or 
calculated with use of an accurate expression for the drag. For rough estimates of the order 
of magnitude of the expected changes in inclination we shall here simply postulate arbitrary 
values of 7,9 and 7,9 and try to guess plausible values of 7; and ry. 


8. A Simplification 


Oblateness is known to produce no secular change in the inclination @, but only an oscilla- 
tion with amplitude less than 0.05°. Suppose we also anticipate a result of the present paper, 
viz, that drag produces only a small secular change in @, of the order 0.3°. For a polar orbit, 
with @=90°, we can then neglect the term in (23.1) that involves cos @. 

For any other orbit this term, although nonvanishing, is still small in comparison with 
the term fm~'p. ‘To show this, note that the ratio ? of the two terms satisfies 


R< mw,r* p. (31) 


With use of (26), it follows that 


If we neglect this term we need consider only the simplified system 


Pp tm"p (23. 1a) 


pe for? cos? y sin @. (23.2) 
9. Stage (2): Quasi-Steady Spiral Orbit 


It is convenient to treat stage (2) first. In this stage we shall assume that the change in 
rin one revolution is so small that we can neglect the corresponding change in f compared to 
that in y. Presumably / is proportional to the density, which in turn presumably varies like 
K exp[—(r—r,)/s], where r,=mean earth radius=radius of sphere of same volume as the earth 

3,959 miles and where s is of the order 10 miles. Our neglect of the variation in f in one 
revolution is thus equivalent to assuming that the corresponding decrease in 7 is somewhat less 
than one mile. (Clearly it would not do to make such an assumption for stage (3), when 7 
may be diminishing much more rapidly.) 

To obtain the secular change, we can then “secularize”’ the right side of (23.2) by replacing 
cos’ ¥ by 's, thereby obtaining, 


pe 5 fwr* sin @. (3.33) 
Division of (33) by (23.la) then gives 


d6/dp=tmw,r*p-* sin @. (34) 
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Note that the unknown drag function has now dropped out from the theory in any explicit 


form. In stage (2) 
" ry, 3D) 
closely enough, so that 
p=mrr, Hb) 
by (26) and (35). Kquations (30) and (36) then vield 
p ni | GM)?’ . 3a) 
For convenience, we now define g, by 
p=GM=4.7.’, 38 
where ¢,=3959 miles, the above mean radius of the earth, and where g, is a corresponding 


effective value of the acceleration of gravity at sea level. free from centrifueal effects. Here 


), = 32.224 ft/sec? 39.1 
r,/ Qe) ?=805.4 sec. 39.2 
Then 
p=Mr. ger. 1() 
From (34) and (40) it follows for stage (2) that 
ese 6 dé—0.25w,7 q r?dr }] 
and thus that 
“A, tan(@,/2) w "7 » 
“ese 6dé—In ~ (r,/q,)°? | (re/r ry) }? 
0 tan(é@,/2 6 ; : 


where 7, and @, are the radius and the inclination of the orbit at the beginning of stage (2) and 


r, and @ are the corresponding quantities at the end of stage (2). 


AN EXampLe: Suppose that the initial altitude in stage (2) is 200 miles and the final altitude 


100 miles. Then 


} 3959 + 200 —4159 miles 433.1 
ro = 3959+ 100—4059 miles. $3.2 
We find 
rol (ry/7 O.03S60. $3.3 
Then since w, = 22/86.164 see~'. it follows from (42). (43.3). and (39.2) that 
tan (@,/2) = 
In 0.000378 1. 44 
tan (6,/2) 


Thus @ does not differ much from 4, If we now put 
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then, by (42) 
tan n»=exp (—Q) tan no. (46) 
To find from (46) a power series in Q for 
Ad =0.—?; (47) 
that will be valid for all 6), 6) in the interval (0, 7), we use the fact that » and are almost equal 
to each other, with each lving in the first quadrant, and proceed as follows. Multiply (46) 


by 7 cos nm and add cos n, to each side. 
The result is 


COS Ny SCC N EXP in=—COS n+? (SIN no) Exp (—Q) (48.1) 
(exp imo)|1—7y sin np exp (—in)], (48.2) 

where 
¥=1—exp (—Q), (48.3) 


a small positive quantity. Then 
n—m = Im In [1—?ysin my exp (—ino)] (49) 
Y SIN Np COS ny) —O.5y7 sin 2p sin? m+... , (50) 
on expansion of the logarithm and selection of the imaginary part. Thus 
Ad =2(n— 1) y sin 6;—y" sin 6, sin? (@;/2)+.... (51) 
With use of (48.3), this becomes 
Ad @ sin 6,+ Q7[0.5—sin® (0,/2)] sin +... , (52) 
which may be shortened, accurately enough, to 


: 1... a e 
A@ ( sin 6, —_ Je) *| (1/1) 2 — ("2/r,) 2] sn 8). (53) 
) 
Then, since 
@, =e . _ Fs 
6 (7,/d,) 0.009789 radians =0.5609°, (54) 
we have, as our final result for stage (2): 
A@ 0.5609[ (4/7) °° — (rar?) sin 6 degrees, (55) 


where 7,=3959 miles. 
For the above example, where the altitude is assumed to diminish from 200 miles to 100 


) 


miles in stage (2), insertion into (55) of the numerical value 0.03860 from (43.3 gives the result 
A@ 0.0217 sin 6, degrees. (56) 


This is a very small effect. We next consider stage (1), where we can often expect a larger 
effect, but where the calculations are more uncertain. 


10. Stage (1): Shrinking Elliptic Orbit 
In stage (1) the density and thus the function f in (23.1a) and (23.2) are appreciable only 
when the satellite is very close to perigee. In effect d@=0 and dp 0 unless r- r,. If we 
divide (23.2) by (23.1a), the result is 


esc 6 dé Me@,r-p~ COs” yidp. (57) 
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In stage (2) we denote initial quantites by the subscript 0 and final quantities by the subscript 


1. The integral of (57) over the whole of stage (1) is then 


tan 0 2 *p ) 9 9 = 
esc @ 14 In me, rp "co yl p. (3) 
f) tan (6,/2 p 


In the light of the above remarks we may put 


r? COs” y rs COS" V> 59) 
for ans appreciable increment dp, with the result 
tan (4,/2 *p 
In mo r2p-? cos? v,dp. (60) 
tan (@,/2 p pl Yet 


Now the perigee radius changes little even in stage (1), so that we shall not be far wrong if 


we replace , in (60 by either yg OY 1). Let us therefore replace ?. by the arithmetic mean 


P=$(r +11). (61) 


In (60) y, is the angle between the line of nodes and the radius to perigee (the line of 


apsides When oblateness is taken into account (but not drag), the rate of change of y, is given 


by Garfinkel (see footnote 2, eq. 53), whose result in the present notation is 
y, Se(r,/1)-n 1.25 sin’ 6 , (62) 
where 


e=(/,—T1,) Mr? =0.00109, (63) 


the quantities 7, and /, being respectively the polar and transverse moments of inertia of 
the earth. We have here replaced Garfinkel’s equatorial radius by our mean radius r,, an 
unimportant change. The quantity / is the semi-latus rectum, @ is the inclination of the 


orbit, and » is the mean angular velocity of the satellite in its orbit: 1. e 


/ Mea ; 64.1 
where the semi-major axis a is given by 

a=Il(1—e (7 Poh (64.2) 

With the use of (64), (62), (38), and (39), we find 
vy, —=3ey,/r )A(r,/a)@(1—e) (1 — 1.25 sin? 6 (65.1) 
20.1 °(r,/a l—« (1—1.25 sin’ @) per day. (65.2) 
As an example, suppose that 7 3959+ 400) miles and 1 (3959-2000) miles and 
let us guess that Pas=9 / (3959+ 200) miles. (This latter figure corresponds to the 
beginning of stage (2), as assumed in our example above.) Then, by (64.2) and (29.2), a9 =5159 
miles and ¢) 0.1552, so that e§=0.0241. Also a,;=4159 miles and ¢,=0. The factor 
ra i—< then increases from 0.415 to 0.842, with a mean value of about 0.628. Thus 
(y, 12.6°(1—1.25 sin’ @) per day. (65.3) 
This amounts to 12.6° per day for an equatorial orbit or —3.15° per day for a polar orbit. 


Thus for a polar orbit the period of ¥, is about 3.7 months. If the lifetime of a polar 


satellite is less than this, it remains doubtful what to do about cos y, in (60), although even 
then we can establish an upper limit for @ of 
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If stage (1) lasts for a vear or more, then ¥, goes through a number of cycles in this state, 
unless sin? @ is too close to the value 0.8, for which the inclination 6=63.4°. If we consider 
various possible orbits and let @ approach 63.4°, the period of ¥, becomes greater and greater, 
until at this critical inclination the angle ¥, would remain constant at its initial value during 
the whole of stage (1). 

For durations of stage (1) longer than a vear and unless @~63.4°, however, we may find 


the secular change in @ by replacing cos* ¥, by 3 in (60), with the result 
tan (6/2) oe —» 1 ' = 
"tan (6/2) ee ae ae — 
= —().5w,(7, Je) AF rr. /ry)?2— (r /bo) 7], (67) 
with the use of (25) and (38). Here 
ls *=0.5(rao' +150) (67.1) 


from (29.1). 


We may evaluate @,—4 from the logarithmic tangent expression, in the same manner 


that we found @,—@, in stage (2). From (67) we find 


A,—86, —0.5 LL (0#63.4°) (68) 


and from (60) we find 


6,—@ /. cosy , (6@+63.4°. with y, effectively constant). (69) 
In (68) and (69) 
L=wo,(r./g)°(/ro* (rer (r,/ly)”*]. (70) 
In anv case 
d,—0,|\< L. (71) 


An Exame.e: We consider the case where the initial perigee and apogee altitudes are respec- 
tively 400 miles and 2000 miles and assume that the final perigee and apogee altitudes of stage 
(1) are both 200 miles. The end of stage (1) in this example is then the same as the beginning 
of stage (2) in the previous example, so that the two examples together constitute one com- 


pound example. The 7,0=— 3959+ 400 =4359 miles, 7.9=—3959+ 2000= 5959 miles, Ypi~Fat 
r, — 3959+ 200=4159 miles, so that L=0.346°. Then, by (68) through (71), we have for 
Stage (1) (In miles, perigee 400->200; apogee 2000-+200) 
A a 0.173° sin @ (0) © 63.47) (72) 
: - a ‘ . ° - 
0,—8, 0.346 cos? y,sin 4, (A, so close to 63.4° that w, Is effectively constant.) (73) 
In any case 
8,—0, <0.346° sin 4. (74) 


These results ma be compared with those for 
Stage (2) (In miles, altitude 200—>100) 
9.—8; 0.022° sin 6; = —0.022° sin % 
The over-all results are that 


6.—86 0.195° sin @ (0) # 63.4°) 


9,—O, <0.368° sin %, In any case. 
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ll. Summary 


The motion of an earth satellite relative to the earth’s atmosphere produces a secular 
change in the inclination of the orbit. In predicting this change it is assumed that atmospheric 
drag acts in the direction of the air velocity relative to the satellite and that its magnitude 
depends on the relative air speed and in such a way on the altitude that it is appreciable only 


near the perigee, in the elliptic stage of motion. No other assumption is made about the drag. 

If @ is the inelination of the orbit, p the orbital angular momentum, m the mass of the 
satellite, w, the angular speed of rotation of the earth, r the distance of the satellite from the 
earth’s center, and py the angle between the radius vector and the line of nodes, it is then 
shown that ese @ d@=mw,r°cos*ydp, to a good approximation, independently of the magnitude 
of the drag. 

The motion is separated into an initial elliptic stage, a quasi-steady spiral stage, and a 
final ballistic stage. The secular change in @ is deduced for the spiral stage by integrating the 
above differential equation with use of an average value for r and the average value 's for 
cos? y. For the elliptic stage the secular change is found by putting r*cos*y equal to the perigee 
value r*,cos* y, in the above equation and then using some results of Garfinkel (see footnote 2) 
concerning the motion of the perigee. Nothing can be said about the final ballistic stage, 
however, without accurate knowledge of the drag. 

The secular change A@ in the inclination, a quantity difficult to measure, is thus expressed, 
for the first two stages of motion, in terms of easily measurable satellite distances. As an 
example, it is assumed that a satellite, with initial perigee and apogee altitudes of 400 miles 
and 2,000 miles, falls to an altitude of 200 miles at the end of the elliptic stage and an altitude 
of 100 miles at the end of the quasi-steady spiral stage. If these data should happen to be 
correct, then 4@ would be —0.022° sin 6 for the spiral stage and —0.173° sin 6) for the elliptie 
stage. The latter figure depends on the assumption that 6 is not close to the critical value 
63.4°, for which there is no motion of the perigee relative to the line of nodes. In any case 
the total A@ would always be less than 0.368°, for the assumed satellite altitudes, before the 
ballistic stage sets in. Onee the ballistic stage begins, no predictions can be made by the 


methods of this paper. 


WASHINGTON, September 15, 1958. 
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